Abstract. IR spectroscopy provides a valuable tool 
for the characterisation and identification of interstellar 
molecular species. Here, we present 6-9 fim spectra of a 
sample of reflection nebulae, Hll regions, YSOs, evolved 
stars and galaxies that show strong unidentified infrared 
bands, obtained with the SWS spectrograph on board 
ISO. The IR emission features in this wavelength region 
show pronounced variations. 1) The 6.2 /im feature shifts 

■ from 6.22 to 6.3 /im and clearly shows profile variations. 
,; 2) The 7.7 /im complex is comprised of at least two sub- 

peaks peaking at 7.6 and one longwards of 7.7 fim. In 
, some cases the main peak can apparently shift up to 8 

fim. Two sources do not exhibit a 7.7 fim complex but 
I ^ , instead show a broad emission feature at 8.22 /im. 3) The 

8.6 /im feature has a symmetric profile in all sources and 
, some sources exhibit this band at slightly longer wave- 

■ lengths. For the 6.2, 7.7 and 8.6 /xm features, the sources 
, have been classified independently based on their profile 

■ and peak position. The classes derived for these features 
^ ,. are directly linked with each other. Sources with a 6.2 

fim feature peaking at 6.22 /im exhibit a 7.7 /im com- 
, plex dominated by the 7.6 /im component. In contrast, 
l/^ ' sources with a 6.2 /im profile peaking longwards of 6.24 
fim show a 7.7 fxm complex with a dominant peak long- 
wards of 7.7 fim and a 8.6 fim feature shifted toward the 
, red. Furthermore, the observed 6-9 /im spectrum depends 
on the type of object. All ISM-like sources and a few PNe 
Oh, and Post-AGB stars belong to the first group while iso- 
Q ■ lated Herbig AeBe stars, a few Post-AGB stars and most 
|H ' PNe belong to the second group. We summarise existing 
^ • laboratory data and theoretical quantum chemical calcu- 
. . I lations of the modes emitting in this wavelength region of 
^ • PAH molecules. We discuss the variations in peak position 
and profile in view of the exact nature of the carrier. We 
5J] attribute the observed 6.2 /im profile and peak position 
5^ I to the combined effect of a PAH family and anharmonic- 
ity with pure PAHs representing the 6.3 fim component 
and substituted/complexed PAHs representing the 6.2 fim 
component. The 7.6 fim component is well reproduced by 
both pure and substituted/complexed PAHs but the 7.8 
/im component remains an enigma. In addition, the exact 
identification of the 8.22 fim feature remains unknown. 
The observed variations in the characteristics of the IR 
emission bands are linked to the local physical conditions. 
Possible formation and evolution processes that may in- 
fluence the interstellar PAH class are highlighted. 
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1. Introduction 



tema et al. 1996 ; Molster ct al. 1996 ; Roelfsema et al 
1996| ; [Moutou ct al. 19994 H; |Peeters et al. 1999]) . It was 



Mid-inuaiud spetiLia uf many suuiues aiu dumiuatud by 
the well-known emission features at 3.3, 6.2, 7.7 and 
11.2 /im, commonly called the unidentified infrared (UIR) 
bands |( cf. i GlUeLL cL al. 1973|, i Gcballe cL al 



recognised some time ago that the 7.7 /im complex ap 
pears either with a dominant 7.6 fim component or with 
the dominant component peaking at 7.8-8 fiui (Bregman 



et al. 1 986h. Tliusc UIR bands aie associated with a wide 



I Culnjii uuiiiiiiciiiL (juiiipuiieiiL peuis-iii^ cll t .o~o pin [ ^lk^^lllchl. 

±joo^ | v^uiit.ii | -^Qgg , Q^Yien et al. 19891 ). In addition, it was found that 



variety of objects - including HII regions, Post-AGB stars, 
PNe, YSOs, the diffuse ISM and galaxies - and are gener- 
ally attributed to polycyclic aromatic hydrocarbon (PAH) 
molecules (Leger &: Puget 1984 ; AUamandola et al. 1985| ; 



Puget fc Leger 198S ; AUamandola et al. 1989| ), although 



the exact molecular identification of the carriers remains 
unknown. Beyond serving as simple PAH indicators, they 
can serve as red-shift indicators, as tracers of elemental 
evolution in external galaxies, as tracers of chemical evo- 
lution and can be used to probe environmental conditions 



within the objects ( 


Genzel et al. 1998; 


Lutz et al. 1998; 


Helou 1999; 


Serabyn 1999^ Genzel & Cesarsky 200C|; 


Helou 


et al. 2000; 


Joblin et al. 2000; 


Hony et al. 2001; 


Vermeij 


et al. 2002; 


Verstraete et al. 2001 


)• 



The region from 6 to 9 fj,m reveals a number of emis- 
sion features with bands at 5.2, 5.7, 6.0, 6.2, 6.8, "7.7" and 
8.6 /zm. The "7.7" /zm feature is particularly important 
as it is the strongest of the interstellar "UIR" bands and, 
as such, can be used to probe objects in which the other 
features are weak. 

Until quite recently, most of the interstellar emission 
bands were considered to be more-or-less invariant in po- 
sition and profile. Although some minor variations were 
noted, by and large the 6.2 /im feature was considered 
fixed at 6.2 /im, regardless of the reported shift in peak 
position by [Molster et al. (1996| ). The "7.7" fxm band 
was generally treated similarly in spite of earlier papers 
showing this band is comprised of at least two variable 
components (e.g. Bregman 1989; Cohen et al. 1989; Bein-| 
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Based on observations with ISO, an ESA proioct with in- 



struments funded by ESA Member States (especially the PI 



countries: France, Germany, the Netherlands and the United 



Kingdom) and with the participation of ISAS and NASA 



the former profile is associated with HII regions and the 
one peaking near 7.8 /im is associated with planetary neb- 
ulae (Bregman 1989; Cohen et al. 1989| ). Recently, thanks 
to the high resolution spectra obtained with ISO, more 
subpeaks of the 7.7 /im complex were reported near 7.2 
to 7.4 and 8.2 fim ( [Moutou et al. 1999a] Jb|). 

In Sect. ^ our sample and the observations are pre- 
sented, the data reduction, the influence of extinction and 
the decomposition of the spectra are discussed. Sections 
^J, ^^2 and 3.3 analyse the 6.2, 7.7 and 8.6 /im features re- 
spectively. The link between the observed variations in the 
6.2, 7.7 and 8.6 /im features and the connection with the 
type of object is highlighted in Sect. IJ. Section || presents 
the observed trends. The spectral characteristics of PAHs 
in this wavelength range as measured in the laboratory 
and calculated by quantum chemical theories arc sum- 
marised in Sect. ||. Section |^ highlights the astronomical 
implications. Finally, in Sect, ^jour main results are sum- 
marised. 



2. The data 

2.1. Sample 

The sample includes 57 sources from a wide variety of ob- 
jects, ranging from Reflection Nebulae (RNe), HII regions, 
Young Stellar Objects (YSOs), Post-AGB stars, Planetary 
Nebulae (PNe) to galaxies (see Table ||) . We give in Table 
characteristics of the sources; i.e. the extinction in the 
K-band, A^, the spectral type of the illuminating source, 
and an estimate of the incident UV flux density at lOOOA, 
Go , at the location where the PAH emission originates in 
units of the average interstellar radiation field {Xu\=1.6 
10-6 W/m2, [Habing 196^). 
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Table 1. Journal of observations. The coordinates of the SWS pointing are given. 
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Source 


a 

(J2000)" 


5 

(J2000)'' 


TDT*" 


Obs. 
mode^ 


Rof. 




Sp. Type'' 


Gg 


Object Typo 


NGC 253 


00 


47 


33 


19 


-25 


17 


17.20 


24701422 


01 


4) 


1 




- 


- 


Scyfcrt Galaxy 


W 3A 022r9+6125 


02 


25 


44 


59 


+62 


06 


11.20 


64600609 


01 


2) 


2 


1.5 


06 


1E4 


CHII 


IRAS 02575+6017 


03 


01 


31 


28 


+60 


29 


13.49 


15200555 


01 


2) 


2 


2 




1E5 


CHII+YSO 


IRAS 03260+3111 


03 


29 


10 


37 


+31 


21 


58.28 


65902719 


01 


3) 


3 




B9 


2E4 


non-isolatcd Herbig Ac Be stars 


Orion PKl 


05 


35 


13 


67 


-05 


22 


08.51 


68701515 


01 


4) 


4 


0.15 


06 




HII 


Orion PK2 


05 


35 


15 


79 


-05 


24 


40.69 


83301701 


01 


4) 


- 




06 




HII 


OrionBar D8 


05 


35 


18 


22 


-05 


24 


39.89 


69501409 


01 


2) 


5 




06 




HII 


OrionBar BRGA 


05 


35 


19 


31 


-05 


24 


59.90 


69502108 


01 


2) 


- 




06 




HII 


OrionBar D5 


05 


35 


19 


81 


-05 


25 


09.98 


83101507 


01 


2) 


- 




06 


5E4 


HII 


OrionBar H2S1 


05 


35 


20 


31 


-05 


25 


19.99 


69501806 


01 


4) 


6 




06 


7E3 


HII 


OrionBar D2 


05 


35 


21 


40 


-05 


25 


40.12 


69502005 


01 


2) 


- 




06 




HII 


NGC 2023 


05 


41 


38 


29 


-02 


16 


32.59 


65602309 


01 


3) 


7 




B1.5V 


3E2 


RN 


HD 44179 


06 


19 


58 


20 


-10 


38 


15.22 


70201801 


01 


4) 


8 




B8V 


5E6 


Post-AGB star 


IRAS 07027-7934 


06 


59 


26 


29 


-79 


38 


48.01 


73501035 


01 


2) 


9 




WCIO 


2E7 


PN 


M 82 


09 


55 


50 


70 


+69 


40 


44.40 


11600319 


01 


4) 


1 




- 


- 


starburst galaxy 


HR 4049 


10 


16 


07 


56 


-28 


59 


31.31 


17100101 


01 


2) 


8,10 




B9.5Ib-II 




Post-AGB star 


IRAS 10589-6034 


11 


00 


59 


78 


-60 


50 


27.10 


26800760 


01 


2) 


2 


1.5 




1E5 


CHII 


HD 97048 


11 


08 


04 


61 


-77 


39 


18.88 


61801318 


01 


4) 


11 


0.12 


AO 


1.7E4 


non-isolatcd Herbig Ae Be star 


HD 100546 


11 


33 


25 


51 


-70 


11 


41.78 


27601036 


01 


1) 


12 


0.03 


B9Vne 


9E3 


isolated Herbig Ae Be star 


IRAS 12063-6259 


12 


09 


01 


15 


-63 


15 


54.68 


25901414 


01 


2) 


2 


1.5 




1E5 


CHII 


IRAS 12073-6233 


12 


10 


00 


32 


-62 


49 


56.50 


25901572 


01 


2) 


2 


1.5 


06-07.5 


1E6 


CHII/star forming region 


IRAS 13416-6243 


13 


46 


07 


61 


-62 


58 


19.98 


62803904 


01 


3) 


- 








Post-AGB star 


circinus 


14 


13 


09 


70 


-65 


20 


21.52 


07902231 


01 


4) 


13 




- 


- 


Scyfcrt 2 galaxy 


HE 2-113 


14 


59 


53 


49 


-54 


18 


07.70 


43400768 


01 


2) 


14 




WCIO 


6E4 


PN 


IRAS 15384-5348 


15 


42 


17 


16 


-53 


58 


31.51 


29900661 


01 


2) 


2 


1.5 




5E4 


CHII 


G 327.3-0.5 


15 


53 


05 


89 


-54 


35 


21.08 


11702216 


01 


1) 


- 




38,000 




HII 


IRAS 15502-5302 


15 


54 


05 


99 


-53 


11 


36.38 


27301117 


01 


2) 


2 


3.1 




3E6 


CHII 


IRAS 16279-4757 


16 


31 


38 


20 


-48 


04 


06.38 


64402513 


01 


3) 


15 








Post-AGB star 


CD -42 11721 (off) 


16 


59 


05 


82 


-42 


42 


14.80 


28900461 


01 


2) 


3,16,17 


0.7 


BO 




non-isolatcd Herbig Ac Be star 


CD -42 11721 


16 


59 


06 


79 


-42 


42 


07.99 


64701904 


01 


2) 


16,17 


0.4-0.7 


BO 




non-isolatcd Herbig Ac Be star 


IRAS 17047-5650 


17 


00 


00 


91 


-56 


54 


47.20 


13602083 


01 


3) 


9 




WCIO 


5E6 


PN 


IRAS 16594-4656 


17 


03 


09 


67 


-47 


00 


47.90 


45800441 


01 


1) 


18 




B7 




Post-AGB star 


IRAS 17279-3350 


17 


31 


17 


96 


-33 


52 


49.30 


32200877 


01 


2) 


2 


2.2 




5E3 


CHII 


IRAS 17347-3139 


17 


36 


00 


61 


-31 


40 


54.19 


87000939 


01 


3) 


19 






8E5 


PN 


XX-OPH 


17 


43 


56 


42 


-06 


16 


08.00 


46000601 


01 


4) 


- 




Ape 




variable star, irregular type 


Hb 5 


17 


47 


56 


11 


-29 


59 


39.70 


49400104 


01 


3) 


17 




120,000 




PN 


IRAS 18032-2032 


18 


06 


13 


93 


-20 


31 


43.28 


51500478 


01 


2) 


2 


1.1 




2E5 


CHII 


IRAS 18116-1646 


18 


14 


35 


29 


-16 


45 


20.99 


70300302 


06 




2 






8E4 


CHII 


GGD -27 ILL 


18 


19 


12 


03 


-20 


47 


30.59 


14900323 


01 


2) 




2 


Bl 


1E6 


star forming region 




18 


19 


12 


00 


-20 


47 


31.10 


14802136 


01 


2) 


2,20 










MWC 922 


18 


21 


16 


00 


-13 


01 


30.00 


70301807 


01 


2) 


- 




Be 


6E6 


emission-line star 


IRAS 18317-0757 


18 


34 


24 


94 


-07 


54 


47.92 


47801040 


01 


2) 


2 


2.0 


08 


1E5 


CHII 


IRAS 18434-0242 


18 


46 


04 


09 


-02 


39 


20.02 


51300704 


06 




21 


1.6 


03-05 


2E6 


CHII 


IRAS 18502+0051 


18 


52 


50 


21 


+00 


55 


27.59 


15201645 


01 


2) 


2 




07 


1E6 


CHII 


HD 179218 


19 


11 


11 


16 


+ 15 


47 


18.58 


32301321 


01 


3) 


22 


0.37 


B9 


>2E4 


isolated Herbig Ac Be star 


IRAS 18576+0341 


19 


00 


10 


50 


+03 


45 


47.99 


32401203 


01 


1) 


23 




15,000 


4E3 


LBV 


BD +30 3639 


19 


34 


45 


19 


+30 


30 


58.79 


86500540 


01 


3) 


14 




WC9 


1E5 


PN 


IRAS 19442+2427 


19 


46 


20 


09 


+24 


35 


29.40 


15000444 


01 


2) 


2,20 




07 


7E6 


CHII 


BD +40 4124 


20 


20 


28 


31 


+41 


21 


51.41 


35500693 


01 


3) 


24 


0.3 


B2V 


1E4 


non-isolatcd Herbig Ae Be star 


S 106 (IRS4) 


20 


27 


26 


68 


+37 


22 


47.89 


33504295 


01 


2) 


25 


1.4 


08 


2E5 


YSO 


NGC 7023 I 


21 


01 


31 


90 


+68 


10 


22.12 


20700801 


01 


4) 


7 




B3 


5E2 


RN 


CRL 2688 


21 


02 


18 


79 


+36 


41 


37.79 


35102563 


01 


3) 






F5Iae 


5E3 


Post-AGB star 


NGC 7027 


21 


07 


01 


70 


+42 


14 


09.10 


55800537 


01 


4) 


8 




200,000 


2E5 


PN 


IRAS 21190+5140 


21 


20 


44 


89 


+51 


53 


26.99 


74501203 


06 




21 


0.0 




6E5 


CHII 


IRAS 21282+5050 


21 


29 


58 


42 


+51 


03 


59.80 


05602477 


01 


2) 


8 




09 


1E5 


Post-AGB star 


IRAS 22308+5812 


22 


32 


45 


95 


+58 


28 


21.00 


17701258 


01 


2) 


2,20 




07.5 


1E3 


CHII 


IRAS 23030+5958 


23 


05 


10 


60 


+60 


14 


40.99 


75101204 


06 






0.0 


06.5 


8E3 


CHII 


IRAS 23133+6050 


23 


15 


31 


39 


+61 


07 


08.00 


56801906 


01(2) 


2 


0.4 


09.5 


3E5 


CHII 



Units of Q ar e hours, minutes, and seconds, and units of S are deRre e s, arc minut e s, and arc seconds. : each ISO o bservation 



is given a unique TDT (Target Dedicated Time) number. SWS observing mode used (see ie Graauw et al. 1996 ). Numbers 
in brackets corr espond to the scann ing s peed. see text for details. 

References : 1 : ISturm et al. (2000|); 2 : [Peeters et al. (2002|); 3 : |Van Kerckhoven et al. (2000|); 4 : [Rosenthal et al. (19991); 5 
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Beintema et al. 



Szczerba et al 



Cesarsky et al. (20001); 6 : IVerstraete ct al. (20011); 7 : iMoutou et al. (1999a|); 8 
(2001 ^10 : iMolster et al. (19961); 11 : |Van Kerckhoven et al. (1999|); 12 : [Waelkens et al. (19961); 13 : Moorwood et al. (W9(iy, 



14 : [Waters et al. ( 19981) ; 15 : [Tielens et al. ( 199E|) ; 16 : [Benedettini et al. ( 199 j) ; 17 : |Hony et al. ( 200l|)Tl8 : parcfa-Lario ct al 



(1999^ 
et al. ( 



19 



(1996[) 



21 



Cohen et al. (1999[); 20 : [Roelfsema et al 
00C|); 24 : [Wesselius et al. (199e| ); 25 : [van den Ancker et al. (2000[ ) 



Peeters et al. (1999); 22 : Meeus et al. (2001); 23 : Hrivnak 



4 



E. Peeters: The rich 6 to 9 fim spectrum of intersteUar PAHs 



For the compact HII regions (CHII) present in our et al. 1996). The spectra were taken using the AOT 01 



sample, Martin- Hernandez et aL (2002a) estimated Ak 
based upon HI recombination hnes to be between and 
2.7 magn. Ak is tak en from [Cidale et aL (2001 ) for 
CD -42 11721, fr om Miroshnichenko et al. (1999|) for 



HD 1792 18, from [Everett et aL (19951) for Orion Peakl 



and from |van den Ancker (19991) for IRAS 03260, GGD - 
27, S 106, HD 97048, BD +40 4124 and HD 100546. 

For most sources, the spectral types are taken from 
Simbad. The s pectral type of IRAS 12073 and IRAS 18434 
are taken from Kaper et al. (2002a| ,|b|, and private commu- 



nication) and that of IRAS 16594-4656 is from 3u et al. 
(20010 The effective temperatures for Hb 5, NGC7027, 



IRAS 18576 and G327 are taken from pesicki fc Zijlstra 
(2000 |)] [Latter et aL (2000| ), |Ueta et al. (200 1| ) and [Ehrcn- 
freund et al. (1997D resp"ectively. 

For the CHII regions and GGD -27 ILL, we have de- 
rived Go values from the observed IR flux and the angular 
size of the PAH emission region (cf. Hony et al. 2001). 
This estimate is based on the assumption that all the UV 
light is absorbed in a spherical shell with the angular di- 
ameter of the HII region and re-emitted in the IR. We 
have used for the size of the HII regions the measured ra- 
dio sizes. This is reasonable since the PAHs are expected 
to be destroyed inside the HII region. The IR flux was de- 



rived from the L/jj given by Peeters et al. (2002) and the 
radio sizes used are taken from peeters et al. (2002 ) and 



Martfn-Hernandez et al (2002b|). The Gp values are sim- 



ilar to those derived by Hony et al. (2001) for the sources 
present in both samples. For the Orion bar, we refer to 
Tielens et al. (199^ ) and Joblin et al. (199^ ) for the given 
Go values. We have taken Gq values for the Herbig Ae Be 
stars from Van Kerckhoven et al. (2002, in preparation) 
who derived Gq from the UV flux between 6 and 13.6 
eV, Fjjv, and the spatial distribution of the PAHs in the 
sources. Fjjv is derived from the observed stellar flux and 
the known spectral type. CRL 2688 has an effective tem- 
perature of ~6400K. Hence, the FUV luminosity is 0.04% 
of the total luminosity of the sta r. The star's lumi nosity 
and the NIR size are taken from [Goto et al. (2002| ). The 
FIR flux of IRAS 17347 and IRAS 18576 are obtained 
by integrating the modified blackbody that is fitted to the 
SWS spectr a. The size of IRAS 17 347 a nd IRAS 18576 are 
taken from pVIeixner et al. (1999| ) and |Ueta et al. (2001]) 
respectively. For MW C 922, the diameter is taken from 
Meixner et al. (1999) and its FIR flux is derived by inte- 
grating the combined SWS and LWS spectrum longwards 
of 20 /im. For the RNe, PNe and Post-AGB stars not men- 
tioned in this paragraph, the Gq values are taken from 
Hony et al. (2001]). 



2.2. Observations 

All spectra presented here were obtained with the Short 
Wavelength Spectrometer (SWS, de Graauw et al. 1996 ) 
on board the Infrared Space Observatory (ISO, Kessler 



scanning mode at various speeds or the AOT 06 mode, 
with resolving power (A/AA) ranging from 500 to 1600. 
See Table |l| for details of the observations. 

2.3. Reduction 

The data were pr ocessed with the SWS Interactive Analy- 
sis package lA^^ ( |de Graauw et al. 1996 ) using calibration 
files and procedures equivalent with pipeline version 7.0 
or later. Further data processing consisted of bad data 
removal and rebinning with a constant resolution. The 
(sub-)features discussed here are present in all available 
scans. 
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1650 1600 1550 1500 
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Fig. 1. 



6.2 6.4 6.6 

Wavelength (/U,m) 
Both the - independently reduced - up and down 

scans of HD 44179 are shown with their respective con- 

tinua in panel a. Panel b shows the normalised profiles of 

the up and down scan. 



In case of high fiuxes, the obtained spectra can suffer 
from memory effects. These memory effects can infiuence 
the general shape of the continuum as well as the profile 
of broad features. The sources in our sample for which 
memory effects are present, are indicated in Table ^ At 
the time the data reduction was done, no memory correc- 
tion tool was available. Hence, in case of memory effects, 
the average of the up and down scans is taken. In order to 
investigate the influence of memory effects on this study, 
we analyse the source that suffers the most from memory 
effects in our sample, i.e. HD 44179, by comparing the up 
and down scans in the region of interest, i.e. 5.5-9 ^m. 
Fig. |l| shows the influence on the 6.2 /im profile. The dif- 
ferences are small, even in this most extreme case. The 
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Two examples to show the richness of the 6-9 fiui 
We recognise the 6.2 and 8.6 features. The 



Fig. 3. 

region. 

" 7.7" H- n feature breaks up in three components; two clear 
bands £ t 7.6 and 7.8 /im and a shoulder at 7.4 /im. Fur- 
thermo:'e, weak features are present at 6.6 and 8.2 /im. 



influence is more severe for the 7.7 /im complex (Fig. 
The blue wing of the feature is affected, as well as the 
relative strength of the 7.6 component. However, the er- 
ror due to detector memory effects (< 5 %) is less than 
the uncertainty on the integrated band intensity. It will 
not hamper the spectral analysis and source classification 
performed in this paper. Hence, it will not hamper the 
analysis done in this paper. Recently, a memory correc- 
tion tool has become available (OLPIO) and, as a check, 
the 



sources suncrmg irom memory cnects nave oeen re- 
reduced. We found that that memory effects do not alter 
significantly the band profiles. In order to be consistent 
with the analysis of the other sources, we did not apply 
this memory correction. 

Two sources in this sample (Orion peak 1 and Orion 
peak 2) have strong atomic emission lines perched on top 
of the 6.2 /im PAH feature. These lines and the PAH 
feature are easily separated at the resolution of the SWS 
instrument. The contribution from any line is removed 
prior to the analysis of the profiles. 

2.4- The spectra 

Fig. 3 shows spectra of two typical sources to illustrate the 
spectral detail present. The complete 6-9 /im spectrum 
reveals an extremely rich collection of emission features 
with bands at 6.0, 6.2, 6.6, 7.0, "7.7", 8.3 and 8.6 nm 



et al. 1999; Verstraete et al. 2001 ). In particular, upon 
close inspection, some of these features are perched on top 
of an emission plateau of variable strength. The beginning 
of this emission plateau seems to be variable and falls 
longwards of 6 /im while it extends until ^ 9 fim. 

From the richness of the region, it is clear that several 
components are present. The well known 7.7 /im feature 
consists of two main features at 7.6 and 7.8 fim plus shoul- 
ders at 7.3-7 .4, 7.45 and 8.2 ^m. For example, Verstraete 



(Beintema et al. 1996; 


Molstcr et al. 1996; 


Roclfscma et al. 


1996; 


Verstraete et al. 1996; 


Moutou et al. 1999a 




Peeters 



et al. (2001 ) fit the total region with several Lorentzian 
profiles and Van Kerckhoven et al. (2002, in prep.) fit the 
7.7 fim complex with 4 Gaussians peaking at 7.5, 7.6, 7.8 
and 8.0 /im. 

2. 5. Decomposition of the spectra 
2.5.1. The continuum 

The profile of the 6.2 fim feature in the spectra of all 
sources is determined by subtracting a local spline con- 
tinuum or a polynomial of order 1 . To assess the sensitiv- 
ity of the resulting profiles to the continuum choice, two 
extreme baselines have been defined and subtracted. In 
general, the influence of the continuum determination on 
the profile is very small and hence does not change sig- 
nificantly the band profiles nor the source classification 
performed hereafter. In some sources however, the con- 
tinuum determination is subject to some freedom. These 
sources are indicated in Table 0. 
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Fig. 4. Illustrative examples of the continua underneath 
the 7.7 and 8.6 /xm features. The dashed line represents the 
general continuum, the full line the second (local) contin- 
uum. See text for details. 



The continuum determination around the 7.7 and 8.6 
/zm features is quite arbitrary. We choose to draw first a 
general continuum splined through points from 5-6 and 9- 
10 /xm and through points near 7 /im, excluding possible 
small features in those regions (see Fig. ^, dashed line) . In 
this way, the influence of a silicate absorption feature in 
some sources (see Table ||) is completely ignored. In addi- 
tion, to separate and study the individual 7.7 and 8.6 /im 
contributions, we have also drawn a continuum under the 
7.7 and 8.6 /im features themselves. This second (local) 
continuum is determined by taking additional continuum 
points near 8.3 /im - between the 7.7 and 8.6 features 
(see Fig. ^ full line). In this way, an underlying plateau 
component is defined. 

Other ways of decomposing the broad, blended bands 
and determining the underlying continuum will yield other 
results. In particular, for different bandshapes (Gaussian, 
Lorentzian, etc.), different continua and profile parameters 



(centr al wavelength and FWH M) are obtained ( Boulanger 
et al. l |j98| ; Uchida et al. 2000 ). However, these differences 



will affect all sources in a systematic way and while this 
will influence the profiles of the derived features, this will 
not affect the source-to-source variations we follow. 



2.5.2. Extinction 

Extinction can have a serious effect on the apparent PAH 
spectrum ( ^poon et al. 200^ , see also Fig. ||). In partic- 
ular, with increasing optical depth of the silicate absorp- 
tion feature, the 8.6 /xm feature is decreased tremendously 
(see Fig. |). To asses the influence of extinction on the 
band profiles and their intensities, the "standard" extinc- 
tion law (^raine 1985|; |Mathis 1990|; [Vlartin & Whittet 



1990) is applied to a template PAH spectrum. This tem- 



plate spectrum is obtained by a continuum divided spec- 
trum of a source suffering no extinction. From the result- 
ing spectrum, the band profiles and band intensities were 
derived in the same way as for the sources considered in 
this paper (see Fig. ||). Although the full PAH spectrum 
changes significantly with increasing Ak, the normalised 
6.2, 7.7 and 8.6 /im band profiles derived with the above 
discussed continuum determination are hardly affected by 
the applied extinction, largely because the extinction is 
quite grey over this wavelength region. But, the derived 
intensities and hence their ratios are certainly influenced. 
For this work, no extinction correction has been applied. 

Five sources show water ice absorption at 6.0 /xm (see 
Table ||) and hence the profile of the 6.2 /xm feature can 
be infiuenced. In view of the profile of the water ice band, 
its influence would be expected to be strongest on the blue 
side of the 6.2 /xm feature. However, the ice absorption is 
very small in our sources and the 6.2 /xm band is situated 
in the red non-steep wing of the water band. Therefore, its 
influence on the 6.2 /im band profile (FWHM and peak 
position) is negfigible (see Spoon et al. 2002| ). 



2.5.3. Normalisation of the profiles 

For comparison, the 6.2 /xm profiles shown in Figs. ^, |^, ||, 
|, |l^ and |ll| are scaled in such a way that the integrated 
fiux within the profile is equal to 1. In this normalization 
procedure, the lower limit (6.1 /im) is chosen so as to ex- 
clude the 6.0 /im emission feature, while the upper limit 
is determined by the most extreme end of the feature (i.e. 
6.6 /im). Analogously, the 7.7 /imband profiles shown in 
Figs. |l^, |lj, |l^ and |l6| were normalised so that the to- 
tal fiux from 7.2 to 8.2 /xm equals one. For the sources 
IRAS 17347, IRAS 07027, He2-113 and HD 100546, we 
normalised the spectra so that the total fiux from 7.2 to 
8.35 /xm equals one in order to cover the 7.7 /im complex 
completely. The 8.6 /im band profiles were normalised so 
that the total fiux from 8.2 to 8.9 equals one. 



3. The band profiles 

In this section, we show how the various types of PAH 
emission spectra found in our sample can be rationalised 
into spectral classes which correspond to different band 
profiles. This is done independently for the three main 
features in the 6-9 /im region : the 6.2, 7.7 and 8.6 /im 



E. Peeters: The rich 6 to 9 fim spectrum of intersteUar PAHs 



7 




6 8 10 12 6.2 6.4 6.6 7.5 8. 8.4 8.8 

Wavelength (/J-m) 

Fig. 5. The influence of extinction on a template PAH spectrum. The top left panel shows the template, a PAH 
spectrum on top of a continuum of 1. The "standard" extinction law is applied to this template spectrum for an Aa' of 
1.5 and 3.0 (middle and lower left panels resp.). As a reference, the template spectrum is plotted in grey. In addition, 
the derived continuum is shown for the extincted spectra. The derived profiles are shown in the right panels. As a 
reference, the profiles of the template spectrum are plotted in grey on top of the derived profiles of the extincted 
spectra. 



bands. Note however that there is some variability of the 
band shapes within a single class. 

3.1. The 6.2 fim feature 

In this section, we classify the 6.2 /xm bands present in our 
sample. In addition, a decomposition of the band profile 
into two symmetric components is discussed. 

3.1.1. The profile of the 6.2 /im feature 

The sources show a pronounced 6.2 /zm feature, some- 
times preceded by a weak feature at about 6.0 /zm. The 
emission profiles of the 6.2 ^m feature are distinctly asym- 
metric with a steep blue rise and a red tail (see Fig. 
Although the profile of the 6.2 /zm feature is similar for 
all the sources, when examined in detail significant differ- 
ences become apparent. A definite range in peak positions 
is present in our sample, varying between 6.19 and 6.29 
/im (see Fig. The width also varies. The peak posi- 
tions and FWHM values for all sources are given in Table 
^. Perusing the derived profiles, we recognise three main 
classes, which we will designate by A, B, and C. 



First, the majority of the 6.2 /xm bands peak between 
6.19 and 6.23 /im. This group will be referred to as class 
A. Note that the strength of the red tail relative to the 
peak strength, and hence the FWHM, varies within this 
class (see bottom panel of Fig. ||). Furthermore, the top 
of the profile can be peaked or rounded off (see top panel 
of Fig. |). 

Second, the remaining sources have peak positions that 
range up to 6.29 /im. We define members of class B as 
those having profiles with a peak position between 6.235 
and 6.28 /im (see Table |^ and Fig. ^ and ^). In general, 
the profiles of class B have a larger FWHM compared to 
those of class A. 

Class B can be further subdivided. In particular, the 
peak position of MWC 922 is shifted toward longer wave- 
lengths compared to class A but it has a similar profile. 
Other members in this group (Bl) contain profiles which 
show a more pronounced red tail than MWC 922 (see Fig. 
^). These variations parallel those in class A. Class B2 is 
represented by HD 44179 in Fig. |^. Their profile is clearly 
less red-shaded compared to class Bl (and class A). Note 
that the profile of HD 44179 shows substructure: one can 
distinguish a blue shoulder. As can be seen from Fig. |l 
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Fig. 6. The normalised PAH CC stretching features of 

class A. The spectra shown exemplify the variations in- 
herent in this class. The vertical dotted lines show the 
range in peak positions in class A. 



the first 'component' occurs at the peak position of class 
A. This is the only source where substructure is revealed 
so clearly. In contrast to the previous classes (Bl and A), 
no variations occur in the relative strength of the red tail 
within this group. Class B3 (exemplified by IRAS 17347 
in Fig. 0) clearly lacks emission before 6.23 /xm and has 
a somewhat higher peak-to-red-tail ratio than B2 sources. 
This subdivision of group B is somewhat arbitrary. Ex- 
amining Fig. ^, there seems to be a gradual progression of 
the peak position to the red from bottom to top. 

Third, IRAS 13416 shows a different profile peaking 
at 6.29 /im. We classify objects with this last profile as 
belonging to class C. Although the emission band in this 
source is very weak and noisy, it is clear that the profile 
is more symmetric than in other classes and its FWHM 
is the smallest found in our sample. The '6.2 /im' band 
in CRL 2688 peaks at 6.29 /im and exhibits the same 
profile as IRAS 13416. Hence, this source is classified as a 
member of class C (see Fig. 

We did not detect in our samp le the 6.35 /im band 
reported by Verstraete et al. (199(: ) in the HII region of 
M17-SW. In fact, a closer inspection of this data suggests 
that the 6.35 /im feature is the result of a cosmic ray hit 
on the SWS detectors. Therefore, we do not include the 
6.35 /im band in our analysis. 

Variations are also seen in the strength of the weak 
6.0 /im feature. It seems to be uncorrelated with the 6.2 
jMa band in both strength and peak position (contrast e.g. 
HD 44179 and He 2-113, cf. Figs | and|l|). 
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Fig. 7. The normalised PAH CC stretching features of 
class B. For ease of comparison, class A and C are repre- 
sented by IRAS 18434 and IRAS 13416. As a reference, 
the profile of HD 44179 is plotted on top of each profile 
in grey. The vertical dotted lines show the range in peak 
positions for class B. 
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Fig. 8. The normalised PAH CC stretching features of 
class Bl. These spectra show the variations inherent in 
this grouping. 



3.1.2. Decomposing the 6.2 /xm feature 

Because the peak position seems to gradually and 
smoothly increase, we hypothesise that all profiles of class 
B consist of a combination of two components with ex- 
treme peak position; one profile of class A and one profile 
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A/Bl 
A 
A 
A 
A 
A 
Bl 
A 
B3 
B3 
A 
A 
A 
A 
Bl 
A 
A 
A 
Bl/2 
Bl 
Bl 
A 
A 
A 
A 
C 
A 
A 
A 
A 
A 
A 



i ntensity fraction of comp. 2 (see Sect. 3.1.2 )^ 



6.212±0.010 
6.223±0.017 
6.227±0.008 
6.216±0.010 
6.207±0.006 
6.215±0.006 
6.206±0.017 
6.222±0.016 
6.225±0.006 
6.212±0.011 
6.217±0.016 
6.214±0.010 
6.268±0.007 
6.268±0.006 
6.210±0.001 
6.260±0.009 
6.223±0.010 
6.221±0.010 
6.251±0.022 
6.217±0.011 
6.205±0.005 
6.299±0.015 
6.210±0.014 
6.255±0.007 
6.222±0.005 
6.228±0.025 
6.211±0.008 
6.219±0.017 
6.227±0.029 
6.224±0.008 
6.212±0.014 
6.246±0.012 
6.215±0.005 
6.259±0.011 
6.270±0.027 
6.223±0.029 
6.209±0.016 
6.227±0.007 
6.205±0.009 
6.243±0.006 
6.224±0.006 
6.215±0.003 
6.214±0.022 
6.257±0.032 
6.249±0.027 
6.239±0.008 
6.218±0.008 
6.203±0.011 
6.210±0.014 
6.213±0.013 
6.290±0.015 
6.213±0.002 
6.210±0.001 
6.213±0.005 
6.205±0.010 
6.205±0.014 
>.217±0.00t 



13.4±0.5 
11.1±0.3 
13.3±2.2 
13.2±0.2 
12.5±0.5 
12.4±0.6 
15.2±0.2 
12.2±0.3 
12.4±0.2 
12.2±0.1 
11.7±1.6 
11.0±1.0 
17.6±0.4 
16.0±1.4 
13.3±0.7 
15.2±0.6 
9.4±0.3 
13.8±0.3 
13.1±1.0 
11.5±0.1 
10.7±1.6 
8.7±1.3 
10.0±2.5 
16.7±0.4 
13.6±0.1 
14.6±1.0 
10.2±0.8 
19.0±0.3 
14.3±2.2 
12.5±0.4 
13.0±0.2 
17.1±0.8 
12.2±0.4 
11.0±4.7 
14.5±4.7 
13.0±2.7 
11.8±0.2 
11.3±1.4 
13.1±0.8 
13.3±0.2 
13.3±0.6 
12.1±0.4 
16.0±0.6 
18.6±3.0 
12.2±0.3 
14.0±0.6 
13.6±0.3 
9.8±0.5 
13.9±0.3 
13.4±1.4 
11.1±1.2 
15.6±0.2 
11.6±1.6 
14.6±0.2 
13.4±0.2 
11.0±2.1 
■^.9±Q.5 



position is do minated by the error on the continuum and is <0.015; based upon the Gaussian fit to the 7.7 /im complex 
(see Sect. 3.2); EPAH/cont [6-9 fim] (see Sect. ^ d: feature detected; ^: strong [NeVI] present on top of the 7.6 /im 
feature, and hampers the determination of the position; the peak intensity of t he 7.6 and 7.8 fim features are equal within 
the estimated error; '': sources with an unusual 8.6 /xm feature (see Sect. 3.3); * and sources suffering from memory 
effects at resp. 6 and 7-9 /im; water ice absorption (6 fim) present; ^: silicate absorption (9.7 /im) present; °: free- 
dom in cont. determination at 6 fj,m; redshifted corrected. Redshifts are taken from the NASA/IPAC Extragalactic Database. 
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Fig. 9. The normalised PAH CC stretching features of 

class C. 
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Fig. 10. Panel a shows the two components used for fitting 
the normalised profiles in class B. Panel b and c show the 
normalised profiles of He 2-113 and MWC 922 respectively 
together with the fit and the two components, reflecting 
their respective contributions. The spectra are interpo- 
lated to the same wavelength grid of resolution 500. 



with the other extreme peak position at 6.29 fim, class 
C. For this two component analysis, we have adopted 
IRAS 18434 as component 1 and IRAS 13416 as compo- 
nent 2. In order to fit a combination of these two compo- 
nents, we interpolated all spectra on the same wavelength 
grid of resolution 500. All profiles of class B are in general 
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Fig. 11. The normalised 6.2 ^m profile of IRAS 18434 

(solid line) with the possible fraction of component 2, i.e. 
class C (dashed line). The dotted line represents the Gaus- 
sian profile derived through subtraction. The spectra are 
interpolated to the same wavelength grid of resolution 500. 
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Fig. 12. The peak position of the 6.2 /im band plotted 
against the fraction of component 2 in the profile (see 
text). The symbols are : HII regions -|-, Herbig AeBe stars 
o, RNe/Be/sym □, Post-AGB stars A, PNe x and galax- 
ies *. 



well fit by this procedure (see Fig. |l^. However, we note 
that for B2/3 spectra, the fit starts slightly short- wards of 
the observed profile. The only exceptions are MWC 922 
and IRAS 17347 which show some subtle differences (see 
Fig. I). Longwards of 6.3 /im, MWC 922 is fit well by com- 
ponent 1 and component 2 seems to be absent. However, 
the peak of MWC 922 is shifted towards the red with 
respect to component 1. IRAS 17347 exhibits a similar 
profile as component 2 with slightly more emission on the 
blue wing. However, the profile of IRAS 17347 is shifted 
to the blue compared to component 2. 
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In fact, the component 1 profile we defined above 
(IRAS 18434) may ah'eady include some contribution of 
component 2 (IRAS 13416). To investigate this case, we 
have subtracted a scaled component 2 from the profile of 
IRAS 18434 (the adopted component 1); a symmetrical 
profile remains peaking at 6.21 and with a FWHM of 0.09 
fim (see Fig. |ll|). Applying this method to all sources, the 
same symmetrical profile remains for all sources. Hence, 
it is possible that the two intrinsic profiles are this de- 
rived symmetric component and the intrinsically symmet- 
ric profile of component 2. The fraction of component 2 
varies from source to source (see Table H) covering the 
full range from 20 to 100 %. The systematic uncertainty 
due to the imperfect template spectra can be estimated 
by degrading the resolution of the template and the fit- 
ted sources from 500 to 100 and is found to be 5%. Based 
upon independent analysis of the up and down scans, we 
have also estimate the uncertainty in the fitting proce- 
dure associated with statistical noise. In Table |^ we quote 
the statistical uncertainty, unless it is less than this es- 
timated systematic uncertainty. Obviously, a strong cor- 
relation is present between the fraction of component 2 
and the peak position of the profile (see Fig. [l^) . The two 
sources whose detailed profile deviates the most in this 
procedure (MWC 922 and IRAS 17347), still agree well 
with the observed trend. It is noteworthy that no source 
in our sample shows this "derived" symmetric profile. One 
source, which we excluded because of its strong ice and sil- 
icate absorption features, does exhibit a 6.22 /im profile 
which closely resembles this derived symmetric profile. It 
will be discussed in a forthcoming paper (Peeters et al. 
2002, in prep.). Furthermore, BD+40 4124 docs exhibit 
a symmetric profile but slightly redshifted with respect 
to the profile derived through the present decomposition 
procedure. 

3.2. The 7.7 fim band 

Although the profile of the 7.7 /xm complex seems similar 
for most sources, as with the 6.2 /xm band, significant dif- 
ferences become apparent when this complex is examined 
in detail. Typically, this band shows major subfeatures at 
7.6 and ~ 7.8 fim (cf. Figs. 3 and ^3|) with possible mi- 
nor subfeatures near 7.3 to 7.4, 8.0 and 8.2 ^m. A definite 
range in relative strength of the 7.6 versus 7.8 /im com- 
ponent is present in our sample, going from a dominant 
7.6 /im component toward a dominant component peak- 
ing longwards of 7.7/im. Furthermore, the 7.7 fim complex 
shifts as a whole. In particular, when the 7.6 component 
is not dominant, the peak position of the whole complex 
varies from 7.79 to 7.97 fim. Whether the peak position 
of the minor 7.6 component also varies in the latter case, 
cannot be determined from the present dataset since it is 
then situated in the wing of the dominant component. In 
this paper, we will make a distinction between those 7.7 
fim complexes with peak position shortwards of 7.7 fim, 



referred to as dominated by the so-called 7.6 fim com- 
ponent and those 7.7 fim complexes with peak position 
longwards of 7.7 fim referred to as dominated by the so- 
called 7. 8 fim feature. This does not necessarily imply that 
the apparent shift in peak position is due to a shift of the 
component itself but can also be due to a different rela- 
tive strengths of the various components giving rise to the 
total 7.7 fim complex. The derived peak positions of the 
7.7 fim complex are given in Table ^. 
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Fig. 13. The typical - normalised - profile of the 7.7 fim 

complex of class A' represented by IRAS 23133, shown 

together with the particular cases within this class. As a 

reference, the typical A' profile of IRAS 23133 is plotted 

on top of each profile in grey. 
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Fig. 14. The normalised 7.7 fim complex of class AB'. As 

a reference, the typical profile of class A' (represented by 

IRAS 23133) is plotted on top of each profile in grey. 
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Fig. 15. The normalised 7.7 fim complex of class B'. 
These spectra show the variations inherent in this group- 
ing. As a reference, the bottom profile represents the typ- 
ical profile of class A' and the 7.7 fim profile of HD 44179 
(a member of Class B'), is plotted on top of each profile in 
grey. The vertical dotted lines represent the nominal 7.6, 
7.8 and 8 /im positions. Clear shifts in peak position of 
the 7.7 fiia complex are present. 



Perusing the derived profiles, we recognise four main 
spectral classes, which we will designate by A', B', AB' 
and C This classification is based upon the relative 
strength of the 7.6 and 7.8 /im features. Sources where 
the 7.6 /im feature is dominant, group in class A'; class 
B' contains sources which show a stronger 7.8 /im fea- 
ture; class AB' groups those sources where the 7.6 and 7.8 
have equal strength. Two sources do not show a 7.7 /im 
complex and hence form class C. 

When looking in detail at the normalised profiles of 
class A', it is clear that almost all sources have the 
same profile. Fig. O shows the possible variations within 



this class. Orion Peakl, Orion Peak2, NGC 7023 and 
NGC 2023 show extra emission between 7.2 and 7.4 /im. 
NGC 2023 has also a more pronounced Ipv.e/Ipy.s ratio. 
IRAS 16279 shows extra emission on the red wing com- 
pared to the typical profile of class A', indicating a con- 
tributor near 8.0 /im. MWC 922 is slightly redshifted as 
well. The latter two sources also show a less pronounced 
ratio of the peak intensities of the 7.6 and 7.8 /im features. 

Two sources, HD 97048 and IRAS 21282, have an 
equally strong 7.6 and 7.8 /im feature (see Fig. |l^) and 
form class AB'. IRAS 21282 is also redshifted compared 
to the typical profile of class A', while HD 97048 has extra 
emission at the red wing near 8.0 /im compared to class 
A'. 

Class B' contains a collection of many different pro- 
files (see Fig. |l^) . All of them are redshifted by different 
amounts compared to the typical profile of class A'. The 
so-called 7.8 /im feature moves from 7.72 towards 7.97 
fim. The strength of the 7.6 fim feature shows large vari- 
ations (cf. NGC 7027 versus IRAS 17347). Whether the 

7.6 /im feature behaves in a similar way as the so-called 
7.8 /im feature concerning peak position, is difficult to say 
since it is situated in the wing of the 7.8 /im feature. 

Two sources, CRL 2688 and IRAS 13416, show no 

7.7 /im complex and no 8.6 /im feature (see §3.3) but 



instead exhibit a similar broad emission feature peaking 
at 8.22 /im (see Fig. These two form group C 

A detailed fit to the profile can be obtained by fitting 
the t otal region with, for exa mple, several Lorentzian pro- 
files ( Verstraete et al. 2001 ). Here we want to obtain an 
estimate of the relative flux contributions of the 2 main 
components rather than to fit the profile. In order to quan- 
tify these variations in the relative strength of the 7.6 and 
the so-called 7.8 /im features, two Gaussians are fit to 
the 7.6 and the 7.8 /im features respectively. In a first 
attempt, the position and the width of these Gaussians 
are fixed to Ao(7.6) = 7.58 /im, FWHM(7.6) = 0.28 /im 
and Ao(7.8) = 7.82 /im, FWHM(7.8) = 0.32 /im. For most 
sources, this works quite well (see Fig. |6|). The peak in- 
tensity of the 7.6 /im band for class A' and the sharpness 
of the top of the 7.7 /im complex is however not well re- 
produced in all cases (see top panel of Fig. |l^) . With this 
method, the extra emission present in four sources (Orion 
Peak 1, Orion Peak 2, NGC 2023 and NGC 7023) on the 
blue side of the 7.7 /im complex is ignored (see middle 
panel of Fig. |l^). In the cases of a few shifted 7.7 /im 
complexes, no good fit is obtained (i.e. for HD 100546, 
HD 44179, He 2-113, BD +30, IRAS 17047, IRAS 07027, 
IRAS 17347). For these sources, a new fit is made in which 
the position of the 7.8 /im Gaussian is treated as a free 
parameter. The FWHM of the two Gaussians is kept fixed 
although - for a single source - it might give a better fit 
if the FWHM were a free parameter. However, to make a 
consistent comparison between sources, we kept the band- 
widths constant. The relative strength of the 7.6 versus 
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7.8 /im feature as well as the intensity ratio of the 7.6 to 
the 6.2 /im feature are summarised in Table |^. Class C' 
is excluded since it has no 7.7 fj,m complex. Based upon 
independent analysis of the up and down scans, we esti- 
mate the uncertainty in the fitting procedure associated 
with statistical noise to be generally less then 10%. An 
estimate of the systematic uncertainty on the integrated 
fluxes can be obtained by applying a different fit proce- 
dure and comparing to our present 2-component decom- 
position. For that reason, we fitted the 7.7 /im complex 
with 4 Gaussians peaking at 7.5, 7.6, 7.8 and 8.0 /im (Van 
Kerckhoven et al. 2002, in prep.). The strength of the com- 
bined 7.5 and 7.6 /im Gaussians in this method can be 
compared to that of the 7.6 /im component applied here. 
Likewise, the combined 7.8 and 8.0 /im Gaussians in this 
method can be compared to the 7.8 /im component ap- 
plied here. This 4-component method provide better fits 
to the spectra with obvious 7.5 /im components (cf. Fig. 
[l6| , middle panel). However, in general, the differences be- 
tween the two methods are small. Hence, these intensity 
ratios are affected by systematic and statistical uncertain- 
ties. In Table H, we quote the larger of these two. 
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Fig. 16. An illustrative example of the two Gaussians fit- 
ted to the 7.7 /im complex. 

Sources dominated by the 7.6 /im feature, always ex- 
hibit a 7.8 /im feature. In contrast, sources with a domi- 
nant 7.8 /im feature, do not always have a clear and dis- 
tinctive 7.6 /im component. Especially, for sources where 
the 7.8 /im feature peaks at the longer wavelength end, 
the 7.6 /im feature is extremely weak - if present (e.g. 



IRAS 17347). Sources exhibiting the broad emission fea- 
ture at 8.22 /im do not show a 7.7 /im complex at all (see 
Fig. 0). 

3. 3. The 8. 6 /im feature 

The 8.6 /im profile derived with the chosen continuum is 
clearly symmetric for all sources (see Fig. |l^ and all of 
them have the same FWHM. Analogous to the 6.2 and 
7.7 /im, a definite range in peak positions is present in 
our sample (see Fig. |l^). Perusing the derived profiles, 
we define 3 main spectral classes A", B" and C". Sources 
with peak position ranging from 8.58 to 8.62 /im group 
in class A", sources with peak position longwards of 8.62 
/im group in class B". CRL 2688 does not exhibit a 8.6 
/im feature and forms group C". IRAS 13416 shows a 
weak feature near 8.6 /im (see Fig. [ij). Careful exami- 
nation of the independently reduced up and down scans 
show that this feature is only present in one of the two 
scans. Hence, we classify this source in group C". Table 
^ gives the class and the peak position derived by fitting 
a Gaussian to the profile. This classification relies heav- 
ily on the adopted continuum. When taking the general 
continuum (see Fig. 0, dashed line) instead of the local 
continuum (see Fig. fffuU line) to derive the profiles, the 
peak position of the 8.6 /im feature shifts for all sources 
slightly toward the blue. As a result, the definition of the 
three classes changes slightly. Two sources previously clas- 
sified in class B" now belong to class A" (i.e. NGC 7027, 
HD 100546). However, most sources remain in their dis- 
tinct classes. Since the presence of a silicate absorption 
feature has no major infiuence on the 8.6 /im profile (see 
Fig. Ijand Sect. 2.5.2), the results discussed here are valid 
for all sources. 

MWC 922 shows an exceptionally strong 8.6 /im fea- 
ture. One HII region, IRAS 18434, has a similarly strong 
8.6 /im feature. The latter source has a strong silicate 
absorption band but correction for the silicate absorption 
will make this band even stronger (see Fig. A detailed 
study of these sources will be presented in Hony et al. 
(2002, in prep.) and Peeters et al. (2002, in prep.) respec- 
tively. 



4. The classes 

In the previous sections, an independent study of the 6.2, 
7.7 and 8.6 /im features was made for our sample. For each 
of the features, different classes were determined. Com- 
paring those classes, an interesting finding is made. The 
class classification of the different bands correlate. This is 
illustrated in Fig. |l^. Sources with a 6.2 /im feature be- 
longing to class A, have a 7.7 /im complex peaking at 7.6 
/im (class A') together with a class A" 8.6 /im feature and 
are referred to as class A sources; while for those with a 
class B 6.2 /tm feature, the 7.7 /im complex is dominated 
by the so-called 7.8 /im component (class B') and their 
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8.6 fim feature is shifted toward the red (class B"). The 
latter sources are referred to as class B sources. The two 
sources showing a single 6.3 /im feature (class C) exhibit 
neither a 7.7 /zm complex nor an 8.6 fim feature (class 
C). Instead, both sources show a broad emission feature 
at 8.22 fj,m. The two sources with an equally strong 7.6 
and 7.8 ^m subfeature (i.e. class AB'), IRAS 21282 and 
HD 97048, exhibit a class A 6.2 /im feature and a class A" 
8.6 fim feature. Note, however, that their 8.6 /im feature 
peaks at the extreme end of class A". Possibly, these two 
sources form an intermediate state between the spectrum 
corresponding to class A and the spectrum corresponding 
to class B. Another example of this type is NGC 7027. 
Its 7.8 fim subfeature is slightly stronger in peak strength 
than its 7.6 fim subfeature, although by a small amount. 
However, this source belongs to class A when consider- 
ing the 6.2 fim feature and class B" concerning the 8.6 
fim feature. Possibly, this is another intermediate state. 
IRAS 18576 and MWC 922 exhibit a Bl 6.2 fxm profile 
while their 7.7 and 8.6 /im features both belong to classes 
A' and A", respectively. Fig. ^provides an overview of 
the PAH spectrum for each of the three main categories. 
Although the classification of the 7.7 /im complex was 
based on the dominant subfeature, i.e. the 7.6 or 7.8 /im 
subfeature, it is clear from Fig. |l^ that the whole spec- 
trum corresponding with class B is shifted compared to 
that corresponding to class A. 

The two sources with the 6.3 /im feature (IRAS 13416 
and CRL 2688; class C), show a similar spectrum com- 
posed of a rather weak 6.3 /im feature, a 8.22 /im fea- 
ture, an extremely weak 11.2 /im feature and a 3.3 /im 
band (see Hony et al. 2002, in prep.). The profiles of the 
8.22 /im emission feature are identical, as are the 6.3 /im 
profiles. The intensity ratio le. a/Is. 22 is however not the 
same. Hence, the 6.3 and 8.22 /im features occur in sim- 
ilar conditions (object type, T, Go, Ue, etc.). But, based 
on these two sources, it seems that no tight correlation, 
in terms of intensities, exist between them. The 8.22 /im 
feature in CRL 2688 and IRAS 13416 is remarkably simi- 
lar in peak position and width to the plateau underneath 
the 7.7 and 8.6 /im features subtracted from the observed 
spectra (see Fig.^ and Sect. 2.5.1). This similarity is of 



course hard to prove. Here we do note that spatial studies 
have shown that this plateau is an independent emission 
component (cf. Bregman 1989; pohen et al. 198*9 ). Two 
sources which are not in our sample, IRAS 22272-1-5435 
and IRAS 07134+1004, both show a very broad feature 
at ~ 8 /im, somewhat shortwards of the 8.22 fim feature. 
This 8 /im profile is much broader then the 8.22 /im fea- 
ture and may be the result of fortuitous convolution of 
a 7.7 /im complex peaking at 7.6 fim and the 8.22 /im 
feature. 

Considering the classification from an astronomical 
point of view, this analysis shows that the PAH spec- 
trum in the 6 to 9 fim region, correlates with the type of 
source considered. All HII regions, reflection nebulae and 



all the extragalactic sources in this sample have a class A 
6.2 fim feature, a 7.7 fim complex peaking at ^ 7.6 /im 
(class A') and a class A" 8.6 fim feature. All Herbig AeBe 
stars that are still embedded in their molecular cloud and 
have an HII region associated with them behave like the 
ISM sources. The isolated Herbig AeBe stars, HD 179218 
and HD 100546, belongs to class B. These two sources 
also show crystalline silicates indicating that disk chem- 
istry may well influence the PAH population. The evolved 
stars are spread over the different classes. Since their out- 
flow is/ was the place of birth of the dust, they likely show 
the evolution of the PAH population in the early phases 
of life. However, they do not reflect - at flrst sight - a 
clear link between the PAH spectrum and the source evo- 
lutionary state. The two Post-AGB stars IRAS 13416 and 
CRL 2688, exhibit the 6.3 and 8.22 fim features (class 
C). The two extremely metal-poor sources HD 44179 and 
HR 4049 both exhibit a class B spectrum while the Post- 
AGB stars IRAS 16279 and IRAS 16594 display a class A 
spectrum. IRAS 21282 is similar to the latter two sources 
except that it has an equally strong 7.6 and 7.8 fim sub- 
feature. The PNe show similar variability in their spectra. 
Hb 5 is the only PN in our sample showing a class A spec- 
trum. IRAS 07027, He 2-113, IRAS 17047, IRAS 17347 
and BD +30 3639 all show a class B spectrum. In this 
class, sources with mixed PAH classes are also present : 
IRAS 18576 and NGC 7027. 

The PAH spectrum in the 6 to 9 fim region apparently 
reflects local physical conditions (object type, T, Go, Uo, 
etc.) or the accumulated effect of processing from the for- 
mation sites in the AGB or post-AGB phases to the ISM 
or the influence of disk chemistry. 

A detailed study on the variations in the PAH spec- 
trum of Herbig AeBe stars and the link with the local 
physical conditions will be given in Van Kerckhoven et al. 
(2002, in prep.). Hony et al. (2002, in prep.) concentrates 
on the PAH spectra of evolved stars and thus on the evo- 
lution of the PAH population in the early phases while 
Peeters et al. (2002, in prep.) focuses on the HII regions. 

5. Correlation studies 

Given the close relation between the main features in the 
6 to 9 fim region, it is of interest to investigate whether 
the strength of these features can shed more light on the 
origin of the observed relations and variations. Since the 
absolute intensities are influenced by the intrinsic lumi- 
nosity and distance of the source, we study the variations 
in the relative strength of the PAH bands. The intensity 
of the 7.6 and 7.8 fim features used here are integrated 
intensities of the Gaussian flts, while the intensities of the 
other PAH features are integrated intensities of the profile. 
The obser ved in tensity ratios can be influenced by extinc- 
tion (see § 2.5.2 ). Since extinction strongly influences the 
intensity of the 8.6 fim feature, we do not consider this 
feature here. 
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Fig. 17. An overview of the possible variations of the main PAH features in the 5-9 /im region, i.e. the 6.2, 7.7 and 
8.6 /xm features. The spectra are normahsed so that the peak intensity of the strongest feature in the 8 yum region 
equals one. Shown is (i) the CHII region IRAS 23133 illustrative of class A (i.e. A, A' and A") with a dominant 
7.6 peak (black line), (ii) the Post-AGB star HD 44179 illustrative of class B (i.e. B, B' and B"), peaking at - 7.9 
/im (dark grey line) and (iii) the Post-AGB star IRAS 13416 representing class C (i.e. C, C and C"), with a broad 
emission feature at 8.22 /im (light grey line). A clear shift in the 7.7 /im complex is present between IRAS 23133 and 
HD 44179. 



Table 3. An overview of the classification of the 6.2, 7.7 and 8.6 /im features. The 6.2 and 8.6 /im bands are 
classified by their peak position (Ap) while the 7.7 /im complex is classified by its domin ant com pone nt, i.e. the 7.6 
fiui component and/or the so-called 7.8 /im component. See text for more details (Sects. 3^, 3_^ and 3^). For each 



spectral class, the range in the local radiation field Go and the spectral types within the class is given. 



Class 

6.2 


/im band 


Characteristics 
7.7 /im band 
comp. 


8.6 /im band 


type of object 


range in Go 


Sp. types 


A A 


~ 6.22 


A' 


7.6 


A" ~8.6 


HII, RN, galaxies, non-isolated 
Herbig Ae Be stars, PN : Hb5, 
2 Post-AGB stars : IRAS 16279, 
IRAS 16594 


3E2 to 7E6 


0,B 






AB 


equal 




IRAS 21282, HD97048 


1.7E4 &1E5 


09 &A0 


B B 


6.24-6.28 


B' 


"7.8" 


B" >8.62 


isolated Herbig Ae Be stars, PNe, 

2 Post-AGB stars : HR4049, HD44179 


6E4 to 2E7 


B,A, 
WC9-10 


C C 


~ 6.3 


C 


8.22 


C" none 


2 Post-AGB stars : IRAS 13416, 
CRL 2688 


5E3 


F5 



When checking for connections between intensity ra- 
tios and object type, we do detect a relation but there 
are always some exceptions present. Nevertheless, when 
plotted against object class (classes A and B), these ex- 
ceptions disappear and clear segregations between the two 
classes become apparent. 



For each class, the average fraction of the total PAH 
fiux, emitted in each feature, its standard deviation and 
the total range present in our sample are shown in Fig. 
p^ . For the total PAH fiux, we only considered the 3.3, 
6.2, 7.7 and 11.2 /xm features. The fraction of the total 
PAH fiux emitted in the 7.6 and the so-called 7.8 /im fea- 
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Fig. 18. Shown is the integrated strength of the 7.8 /im feature versus the 7.6 and 6.2 features, normahsed on 
the total flux emitted in the 3.3, 6.2, 7.7 and 11.2 /im features. Class A is represented by •, class i3 by o and sources 
not belonging to these two classes by A. 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

't.b/'t.s 'y.e/'y.s 



Fig. 20. The relation between the intensity ratio Ir.e/h.s and the peak position of the 6.2 /im band and the 7.8 /im 
component respectively. Extinction is unimportant for this intensity ratio. The symbols are : HII regions +, Herbig 
AeBe stars o, RNe/Be/sym □, Post-AGB stars A, PNe x and galaxies *. 



tures clearly depends on the class of the source (see Fig. 
[l8| , left panel and Fig. [l^). In contrast, the fraction of the 
total PAH flux emitted in the 7.7 /im complex does not 
differ significantly for the two classes. Furthermore, the 
fraction of the total PAH flux emitted in the 6.2 and 11.2 
fim feature differ for the classes A and B while a similar 



fraction is emitted in the 3.3 fim feature for both classes 
(see Fig. |l^, right panel and Fig. ^9|). Therefore, it seems 
that the increase in l6.2/IsPA_ff for class B compared to 
class A is opposite to the behaviour of the 11.2 /im fea- 
ture. However, there is no indication of this when plotting 
l6.2/IsPAff versus Iii.2/IsPyiff . In addition, the fraction 
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6.2 7.6 7. 

PAH Band (/xm) 

Fig. 19. An overview of the fraction of the PAH flux emit- 
ted in each feature for class A and class B. The diamonds 
crossed by solid vertical lines represent the average frac- 
tion of the total PAH flux emitted in each PAH band over 
a given class and its standard deviation. The filled boxes 
represents the observed ranges of the fraction of PAH flux 
emitted in each PAH band within a class. Light grey rep- 
resents class A and dark grey class B. For the total PAH 
flux, we considered the 3.3, 6.2, 7.7 and 11.2 /im features. 

of flux emitted in the 3.3 or 11.2 /im feature does not cor- 
relate with the fraction emitted in the 6.2 /im feature. So, 
this difference in l6.2/IsPA_f/ for the two classes cannot be 
directly linked to only one band (the 3.3, 7.7 or 11.2 /im 
feature). Furthermore, sources with an increased emission 
in the 6.2 /im feature (class B) also have an increased 
emission in the 7.8 ^m feature and hence a decreased 
emission in the 7.6 /im component compared to class A 



(see Fig. 18, right panel). 

In Fig. ^0|, the intensity ratio of the 7.6 to 7.8 /im sub- 
features, i.e. I7.6/I7.81 is plotted versus the peak position 
of the 6.2 /im feature and of the 7.8 /im component re- 
spectively. The peak position of the 7.8 /im component is 
determined from the position of this subpeak in the 7.7 
/im complex itself and not from the applied Gaussian fits. 
The error in this position is given in Table ^ for entrees 
in which the 7.7 /im complex peaks longwards of 7.7 /im. 
For sources with a dominant 7.6 /im component, we esti- 
mate the error to be less then 0.05 /im. It is clear that the 
sources reflect a gradual variation. In contrast, no clear 
correlation between the ly.g/Ir.s intensity ratio and the 
FWHM of the 6.2 /im band is present. 

We also checked for correlations between band strength 
ratios and the local radiation fleld Go , however we do not 
detect any correlations. Furthermore, the local radiation 



field of the sample sources do not show differences between 
the two classes. In addition, we checked for correlations 
between band strength ratios and the integrated band- 
to-continuum ratio, EPAH/cont, in the 6-9 /xm region. 
This ratio traces excitation conditions and/or abundance 
variations. The error on this ratio is dominated by the 
continuum determination and is less then 10%. We took 
the continuum described in Sect. ^; hence, the plateau 
contributes to the continuum. This ratio is given in Ta- 
ble ^. Here also, we do not detect any correlation. But, 
sources with the intensity ratio I7.6/I7.8 < 1 do not have 
a high SPAH/cont ratio (< 0.45) whereas sources where 
I7.6/I7.8 > 1 span a range in EPAH/cont from to 1.2. 

6. The Infrared Emission Features and PAHs 
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Fig. 21. The absorption spectrum of a mixture of neutral 
PAHs (a) compared to the spectrum of the same PAHs 
in their positive state (b). This comparison shows that, 
for PAH spectra, ionisation has a much greater influence 
on relative intensities than on peak frequencies, with the 
features in the 6 to 10 /im region substantially enhanced 
with respect to the rest of the spectrum (Figure adapted 
from AUamandola et al. 199£). 



The IR emission features at 3.3, 6.2, 7.7, 8.6, and 
11.3 /im are now generally thought to arise from vibra- 
tionally excited PAHs due to their similarity with the 
spectra of PAHs taken under conditions which match the 
salient characteristics of the interstellar environment. One 
of the early pivotal results of the laboratory and theoret- 
ical studies on PAHs reported over the last decade is the 
r emarkable effect ionisation has on the infrared spectra 
(ISzczepanski fc Vala 1993| ; [Langhoff 1996| ; |Kim et al. 2001 
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Hudgins & AUamandola 1999a, and rcf. therein). While 
PAH characteristic frequencies are only modestly affected 
by ionisation, the influence on intensity is striking - par- 
ticularly in the 5 to 10 /im region (Fig. |2^). The bands 
grow from the smallest features in neutral PAH spectra 
to become the dominant bands in ionised PAH spectra. 
The CC stretching vibrations grow in intensity because, 
upon ionisation, the charge distribution changes signifi- 
cantly with the CC skeletal vibration, creating a strong 
oscillating dipole whereas the oscillating dipole for the 
CH motions and their corresponding bands are reduced. 
Fig. ^ illustrates this effect, which has been found for 
all PAHs measured in the laboratory to date. However, 
this is mostly limited to PAHs with less than 48 C-atoms. 
Larger PAHs are only now being studied in the labora- 
tory. Preliminary laboratory results for PAHs containing 
up to 60 C-atoms also show this behaviour upon ioni- 
sation. Theoretical studies of large species support this 
enhancement to much larger sized PAHs. For example, 
Bauschlicher (2002) has shown that this holds for C96H24. 
Astronomically, the observation that the 6.2 and 7.7 /im 
features, the focus of this paper, are the most intense of 
the interstellar emission band class, is a clear indication 
that ionisation is an intrinsic characteristic of interstellar 
PAHs. The 3.3 fim feature on the other hand is charac- 
teristic for neutral PAHs. Hence, both ionised and neutral 
PAHs contribute to the interstellar IR emission spectrum. 



(D 




5 6 7 8 9 10 

Wavelength (^im) 

Fig. 22. Comparison of a typical ISM infrared emission 
spectrum (a) with a composite absorption spectrum gen- 
erated by co-adding the individual spectra of 11 PAHs 
(b). The interstellar spectrum is that of IRAS 23133. The 
individual spectra were calculated using experimentally 
measured frequencies and intensities and assigning a 30 
cm^^ FWHM Gaussian band profile, consistent with that 



expected from the interstellar emitters (e.g. Hudgins fc 



TheTE have been many cuiiipai isuiis but ween PAH 
spectra and the interstellar spectra over the years. As the 
observational tools have become more sensitive and the 
laboratory techniques more appropriate to the interstellar 
case, the fits have become more revealing about the 
different PAH populations in different regions and this, 
in turn, has yielded further insight into conditions in 
the emitting regions. Fig. |2^ shows a fit to the emis- 
sion from IRAS 23133 measured with ISO by spectra 
now available in the Ames PAH IR spectral database 
(http : //web99 . arc . nasa . gov/~astrochem/pahdata/ 
index.html). Although the fit is striking, thanks to the 
quality of the new ISO spectra, important differences 
become apparent which shed further light on the inter- 
stellar PAH population. The following differences can be 
seen in going from shorter to longer wavelengths. The 
laboratory band near 6.2 /im falls slightly to the red of 
the interstellar feature, laboratory components centred 
near 7.7 //m do not precisely match the interstellar profile 
at this position, the 8.6 /im laboratory band is weaker 
with respect to the other features in the spectrum than is 
the case for IRAS 23133, and the laboratory component 
near 11.2 /im lies to the red of the interstellar feature. 
The differences in t he 8.6 and the 11 .2 /im regions are 
discussed elsewhere ( Hony et al. 2001 , Janssen, Janssens 
and Swansen 2002, in prep.). Here we focus on the 
interstellar features near 6.2 and 7.7 /im. 



AUamandola 1999b|) . 



Feature assignments in the 5 to 10 /im region 

The 5 to 10 /im region encompasses frequencies which 
originate from a variety of PAH molecular vibrations. 
Pure CC stretching motions generally fall between about 
6.1 /im and 6.5 /im, vibrations involving combinations 
of CC stretching and CH in-plane bending modes lie 
slightly longward, between roughly 6.5 /im and 8.5 /im, 
and CH in-plane wagging vibrations give rise to bands 
in the 8.3 /im to 8.9 /im range. While the well-known 
interstellar features at 6.2, 7.7, and 8.6 /im dominate 
this range, there are at least four weak interstellar 
bands in this region as well, centred near 5.2, 5.6, 6.0 
and 6.8 /im. Their correlation with the major features 
indicates that they too originate from the interstellar 
PAH family. The 5.2 and 5.6 /im features most likely 
correspond to combinations and overtones involving 
the CH out-of-plane fundamental vibrations which 
fall between 11 and 13 /im; the 6.0 /im feature likely 
indicates a carbonyl (>C=0) stretch of an oxygenated 
PAH (a quinone); and the 6.8 /im band probably corre- 
sponds to a weak aromatic CC stretching - CH in-plane 
bending combination mode of PAHs as for example in 
the fluoranthenes, or the aliphatic -CH2- or -CH3 defor- 



E. Peeters: The rich 6 to 9 fim spectrum of intersteUar PAHs 



19 



mation in a methyl or ethyl side-group attached to a PAH. 

There are also other plausible interstellar PAH-related 
species that are likely to be important in the emission 
zones and which should be considered. Some examples 
include PAH clusters and PAH complexes with metals 
such as iron ( met alio cenes ) . Furthermore, one can ques- 
tion whether PAHs remains planar as they grow. And, 
if not, how does the 3-dimensional shape influence the 
IR spectrum? It is noteworthy in this regard that simple 
fuUerenes do not have IR characteristics which coincide 
with any prominent st ructure in the interst ellar emission 
spectra discussed here ( Moutou et al. 1999b ). Further, and 
in spite of much effort, there has not been any report of IR 
active transitions in carbon nanotubes, very large curled 
aromatic networks. If such curled aromatic networks (very, 
very large PAHs) do possess infrared activity, it appears 
to be very weak. Thus here we focus on the PAH molecules 
measured in the laboratory or theoretically calculated. 

The shifts in the profiles shown by the observations 
presented above (see Table ^) are much larger than the 
binsize corresponding to the obtained resolution of the 
data. Hence, they are chemical in nature, arise from local 
excitation conditions and are not due to doppler broaden- 
ing and shifting. Therefore, these shifts and profile vari- 
ations provide important new insight into the variations 
of the interstellar PAH populations in the different en- 
vironments. Interpretation of these new spectral aspects 
require probing deeper into the details of PAH spectro- 
scopic properties in this region than heretofore. Here we 
consider this new information as it affects first the 6.2 /im 
region and then the 7.7 /im region. 



6.1. The position of the PAH CC stretching band near 
6.2 /im .■ experiment and theory 

There are several properties which determine the precise 
peak positions of the infrared active bands which corre- 
spond to pure CC stretching vibrations in any given PAH. 
These include molecular size, molecular symmetry, and 
molecular heterogeneity. The roles each of these play in 
determining the position are discussed below. The charge 
of the PAH molecules also s hifts the band po sition; but 
this trend is not systematic ( Bakes et al. 2001 , Bakes E., 
private communication). In addition, dehydrogcnation in- 
fluences the band position ( Pauzat et al. 1997 ). Ho wever, 
based upon the CH out-of-plane bending modes, Sony 
et al. (2001 ) conclude that dehydrogcnation has little in- 
fluence on the observed interstellar PAH spectrum. 



Molecular size At the small end of the PAH size distribu- 
tion (Cio to '^Cso), the peak wavelength of the dominant 
CC stretching feature decreases steadily with molecular 
size (Hudgins and AUamandola, 1999). The data shown 



in Fig. which has been updated to include more PAHs 
than available in the original study, shows this trend. 
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Fig. 23. A plot of the dominant feature position in the 
PAH CC stretching region as a function of molecular size. 
The positions and typical FWHM of the interstellar 6.2 
and 6.3 components discussed in Sect. 3.1 are in- 
dicated by the horizontal shaded bands. Open symbols 
correspond to experimentally measured positions. Solid 
symbols are frequencies computed using B3LYP quan- 



tum calculations (Hudgins ct al. 2001; Bauschlicher 2002 



Bauschlicher, Hudgins, and AUamandola 2002 (in prep.)). 
Theoretical positions typically have uncertainties of 5-10 
wavenumbers, with a few bands showing shifts of 15 to 20 
wavenumbers. 

However, Fig. |2^ shows that the trend does not con- 
tinue ad-infinitum but seems to die out for molecules with 
more than about 30 to 40 carbon atoms, falling slightly 
longwards of 6.3 /im. This behaviour is consistent with 
the limit of graphite which shows an emission mode at 
6.3 fim ( prainc 1984 ). Thus, as the size grows, the in- 



fluence of boundary conditions on these pure CC stretch- 
ing vibrations within the PAH carbon skeleton vanishes 
when the peak wavelength of this mode in graphite is 
reached. Indeed, it is a general chemical rule that within 
any molecule, the further a given bond or chemical sub- 
group is from the site of a modification, the smaller the ef- 
fect exerted by that change on the properties of the specific 
bond or subgroup. This goes for a molecule's fundamen- 
tal vibrational frequencies as well as its chemical prop- 
erties. For the smallest PAHs, addition of even a single 
ring constitutes a significant modification to the carbon 
skeleton and strongly influences its pure CC stretching 
vibrational frequencies. However, as PAH size increases. 



20 



E. Peeters: The rich 6 to 9 fim spectrum of intersteUar PAHs 



an ever-increasing fraction of the molecule's carbon skele- 
ton lies "far" from the site of any particular modification 
and is thus increasingly insensitive to that modification. 
As a result, the corresponding fundamental CC stretching 
modes that arise within that skeleton are progressively less 
perturbed resulting in a levelling off of their frequencies 
above a certain size (^^30 C atoms as estimated from Fig. 

Importantly from the astrophysical perspective, Fig. 
p3| shows that the maximum wavenumber falls some 20 
cm^^ short of 6.2 /im, the peak position of the Type 
A interstellar bands discussed in Sect. 3.L This mis- 



of the circumcoronene cation by removing first one and 
then two of the peripheral rings, the peak position of the 
strongest pure CC vibration does not shift much if at all. It 
falls at 6.36 /im (1572 cm~^) in both circumcoronene and 
circumcircumcoronene with one ring removed and only 
slightly shifts to 6.33 /im (1580 cm~^) in circumcoronene 
with two rings removed. Thus, we conclude that symmetry 
breaking alone cannot account for the difference between 
the interstellar feature and the PAH spectra considered to 
now. 



match between PAH band position and the interstellar 
feature is only worsened by the approximately 10 cm~^ 



excitec 


PAH f 
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). Thus, 



while the correlation between molecular size with band 
position is confirmed by new experimental and theoreti- 
cal data, for the highly symmetric pure PAHs considered 
here, the strongest, infrared active, pure CC stretching 
mode cannot reproduce the position of the 6.2 /xm inter- 
stellar feature and other factors which can slightly shift 
this frequency must be considered. 



Molecular symmetry Another property which can influ- 
ence IR activity is molecular symmetry. In large, sym- 
metric molecules there are many vibrations which cor- 
respond to very weak or infrared forbidden transitions. 
These vibrations are weak or IR inactive in highly sym- 
metric species because the oscillating dipoles arising from 
atomic motions at one position in the molecule are can- 
celled by identical, but oppositely oriented dipoles aris- 
ing from identical, oppositely phased motions occurring 
elsewhere in the molecule. Reducing molecular symme- 
try eliminates some of this oscillating dipole cancella- 
tion and once IR forbidden modes can become IR active. 
This behaviour is expected independent of the means by 
which symmetry is broken. For large, highly symmetric 
molecules such as circumcoronene (C54H^g, symmetry — 
Dg/i) and circum-circumcoronene (CgeH^^, symmetry — 
Dgft) included in Fig. the calculations indicate that 
there are several infrared inactive CC stretching vibra- 
tions that fall between about 6.15 and 6.25 fj.m, (1626 
and 1600 cm^^) positions which overlap the interstellar 
position. Since there is no reason to expect that interstel- 
lar PAHs would be so highly symmetric, we have explored 
the influence of symmetry breaking on the IR spectrum to 
assess the possibility that asymmetric PAHs would show 
strong IR activity at shorter wavelengths, closing the gap 
between the PAH spectra and the 6.2 /im interstellar fea- 
ture. Several different cases were tested and all show only 
marginal effects (Bauschlicher, Hudgins, and AUamandola 
2002, in prep.). For example, by reducing the symmetry 
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Fig. 24. The effect of PAH nitrogenation on the positions 
of the strongest band in the CC stretching region for the 
three PAHs coronene, ovalene, and circumcoronene. The 
right profile indicates the position for the unsubstituted 
PAH, the left for the nitrogenated species. The position 
of nitrogen substitution is indicated by the filled circle in 
the structures. The vertical grey shaded bands indicate 
the position and width of the interstellar type A and B 
bands discussed in Sect. 3.1. 



Molecular heterogeneity Another means by which IR ac- 
tivity can be induced in otherwise weak or forbidden tran- 
sitions is by the introduction of a hetero atom into the 
symmetric molecular structure. This alters the molecule 
in two ways: it lowers the symmetry and modifies the 
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electronic distribution. To investigate this influence on 
PAH spectroscopy, Bauschlicher, Hudgins, and AUaman- 
dola (2002, in prep.) have studied the effects of introducing 
different hetero-atoms into different positions within the 
carbon network of several PAHs. This affects, in particu- 
lar, the position of the 6.2 fim band. A few examples are 
presented in Fig. Here the positions of the strongest 
band in this region are shown for the cations of the pure 
PAHs coronene, ovalene, and circumcoronene and their 
counterparts which have one of the carbon atoms replaced 
by a nitrogen. In all cases nitrogen substitution shifts the 
strongest CC stretching vibration to the 6.2 /im inter- 
stellar band position. Further within a given PAH, the 
shift depends somewhat on the "depth" of the substitu- 
tion in the carbon skeleton. For example, when nitrogen 
is substituted for a carbon atom on the outermost ring of 
the circumcoronene cation, the strongest band in the CC 
stretching region falls at 6.16 /im; when this substitution is 
made on the next innermost ring (shown in Fig. |^) , it falls 
at 6.21 /im; and when on the innermost ring it falls at 6.19 
/im. Bandshift behaviour when more than one nitrogen is 
substituted is still to be determined. For singly substituted 
PAHs, this peak shift has two causes. First, the insertion 
of nitrogen decreases molecular symmetry, increasing IR 
activity in modes which are normally weak because of the 
higher symmetry as described above. Second, the strongly 
electro negative nitrogen atom modifies the fixed charge 
distribution on the carbon network and induces stronger 
IR activity for these vibrations since the magnitude of 
the oscillating dipole increases. This behaviour can be ex- 
pected independent of the means by which the charge dis- 
tribution is fixed. 

Bauschlicher, Hudgins, and AUamandola (2002, in 
prep.) have also considered oxygen and silicon atom sub- 
stitution. Nevertheless, while these too show a similar in- 
duction of IR activity close to 6.2 /xm, when taking chem- 
ical considerations into account, nitrogen substitution is 
still the most attractive candidate. First, nitrogen can be 
incorporated anywhere within the ring structure without 
compromising the aromatic stability of the tt electron net- 
work. Second, oxygen substitution does not have this ad- 
vantage. Oxygen or O"*" cannot form four bonds and so 
oxygen is not found in the middle rings. Hence, oxygen 
would be most stable only at the edge positions, posi- 
tions which would be more reactive and subject to chem- 
ical attack. Third, although silicon can form four bonds, 
the C-Si bond length is larger then the C-C bond length 
and hence the aromatic rings will be disturbed by sili- 
con substitution. Moreover, silicon has a much lower cos- 
mic abundance then does nitrogen. Thus, at this stage in 
our understanding, large PAHs containing some nitrogen 
seem most plausible to account for the 6.2 /im band po- 
sition. As mentioned before, there are other plausible in- 
terstellar PAH-related species in which one might induce 
activity at this position and these must be investigated 
before a firm conclusion can be drawn. Some examples 



include large PAHs with uneven and irregular structure, 
PAH clusters and PAH complexes with metals such as iron 
(metallocenes). 

6.2. The position of the PAH CC stretching/CH in-plane 
bending feature near 7. 7 : experiment and theory 

The CC stretching/CH in-plane bending vibrations of 
most singly ionised PAHs possess at least one very strong 
feature between about 7.2 and 8.2 //m. Fig. ^ shows a 
graph of the frequencies for the strongest of these modes 
plotted versus the carbon number for the same molecules 
considered in Fig. In this case (Fig. |2^) there is no clear 
relationship between frequency and size in contrast to the 
behaviour for the pure CC stretch (Fig. While a few of 
the smallest molecules do have vibrations which fall at the 
lowest frequencies measured for these modes - suggesting 
such a correlation - similarly sized molecules possess some 
of the highest vibrational frequencies. Furthermore, one of 
the largest PAHs in the laboratory database, hexabenzo- 
coronene (C42H24), has the lowest frequency determined 
to date. Here peak positions for small and large PAHs are 
intermingled. 

In addition. Fig. ^ reveals that there is also a gap be- 
tween about 7.6 and 7.8 fj,m (1316 and 1282 cm~^), with 
a clustering of the data between 1320 and 1355 cm^^ (7.4 
and 7.6 fim). The data in Fig. ^ suggests that the mean 
PAH CC stretching/CH in-plane bending frequency lies 
near 1335 cm~^ (7.5 /im). This behaviour is important in 
view of the results presented in Sect. 3^ which show that 
two components dominate the interstellar emission in this 
region, one peaking near 7.6 /im (class A') , the other long- 
wards of 7.8 nm (class B'). When one takes the roughly 
10 cm^^ red shift into account for emission, it is apparent 
that this family of PAHs can readily reproduce the 7.6 /im 
component, but not the dominant 7.8 /im component. 

As discussed above for the pure CC stretching fea- 
tures near 6.25 /im, we have theoretically considered how 
molecular symmetry and heterogeneity influence the peak 
position of the 7.7 /im complex in order to understand 
what molecular properties might be responsible for the 
interstellar position. Breaking the high symmetry of PAH 
molecules, either by removing rings or by substituting a 
nitrogen atom at several different positions within the car- 
bon network again produces IR activity in modes which 
were weak or completely forbidden (Bauschlicher, Hud- 
gins, and AUamandola 2002, in prep.). In this case how- 
ever, while a few examples are found in which intense 
bands fall in this gap, these results seem random and no 
clear structural relationship has yet emerged. Thus, this 
sort of molecular modification does not account for the 
interstellar 7.8 /im component. 

Furthermore, as with the 6.2 /im feature, dehydro- 
genation and the charge of the PAH molecules do not play 
a m ajor role in determining the peak position (see Sect. 

~T1). 
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Fig. 25. A plot of the dominant feature position in the 
wavelength region of the 7.7 fim complex as a function of 
molecular size. The positions and typical FWHM of the in- 
terstellar 7.6 and 7.8 /im components of class A discussed 



in Sect. 3.2 are indicated by the horizontal shaded bands. 



The two components of the 7.7 /im complex in Class B 
peak at 7.6 and from 7.8 to 8.0 /im respectively. Open 
symbols correspond to experimentally measured positions. 
Solid symbols are fre quencies computed using B3LYP 
quantum calculations (|ludgins et al. 200l|; |Bauschlicher 



2002|7^auschlicher , Hudgins, and AUamandola 2002 (in 
prep.)). Theoretical positions typically have uncertainties 
of 5-10 wavenumbers, with a few bands showing shifts of 
15 to 20 wavenumbers. 



In summary, the 7.6 fim (class A') interstellar band 
is consistent with a mixture of large and small, pure and 
hetero- atomic PAHs. These PAHs also contribute emission 
at 7.8 /im as illustrated in Fig. ^ However, the PAHs so 
far considered reveal a dearth of strong spectral features 
longwards of 7.7 /im and thus have difficulty reproducing 
the class B' 7.7 /im band. There are only a few PAHs in 
the database which have their strongest band at 7.8 /im 
and the data is insufficient to make any generalisations 
concerning the properties which cause strong IR activity 
at this position. As above, there are several other PAH 
related species that are likely to be important in the emis- 
sion zones and which should be considered. 

7. Astronomical implications 

From the wealth of IR spectra, it is clear that the UIR 
bands at 3.3, 6.2, "7.7", 8.6 and 11.2 /im represent a sin- 
gle class of spectral features that come and go t ogether. 
The ISO observations presented here in Sect. 3.1, B.2 and 



p.3| , show that the peak positions and profiles of the 6.2, 
"7.7" and 8.6 /im features vary significantly from source 
to source. Moreover, as discussed in these variations in 
the different bands are correlated with each other and with 
the type of object considered (i.e. HII region, YSO, Post- 
AGB star and so on). These variations stand in marked 
contrast to the behaviour of 3.3 and 11.2 /im bands. While 
the latter modes show some variation in profile, their peak 
position is relat ively invariant with bo t h bands wandering 



by only 0.2% ([Tokunaga et al. 1991^ |Hony et al. 2001 



Peeters et al. 2002, in prep.; van Diedenhoven et al. 2002, 
in prep.). Furthermore, while the 3.3 and 11.2 /im band in - 
tensities are correlated with each other ( Hony et al. 2001 ), 
their strengths do not correlate well with that of the 6.2 
and 7.7 /im features. 

7.1. The observed trends 



It is shown that class A and B have distinctly different 
profiles, peak position and relative intensities of the fea- 
tures in the 6-9/im region. We attribute these differences 
between the two classes to variation in the PAH families 
present around the sources; i.e. different molecular mix- 
tures are found around the sources of the two classes whose 
combined IR spectral fingerprints change the overall PAH 
spectrum in the 6-9 /im region significantly. 

Class B corresponds to regions where Go > 10^. For 
such high Gq, very small grains (VSG) emit in the 6-9 
/im region ([Cesarsky et al. 200C ). Since the range in Go 
present in class B is also present in class A, these VSG 
are not responsible for the observed differences. Further- 
more, these differences in the observed emission features 
and hence in the PAH families is not directly related with 
Go. 



Hony et al. (2001) derived the molecular structure of 
the PAHs based upon the 10-15 fj,m spectra and con- 
cluded that the PNe contain large, compact PAHs with 
long straight edges while HII regions contain smaller or 
more irregular PAHs. Clearly, the HII regions and non- 
isolated Herbig Ae Be stars all belong to class A and have 
a weak I11.2/I12.7 ratio. But, for the PNe considered in 
Hony et al. (2001 ), no link is found between the I11.2/I12.7 
intensity ratio and the profile classes. This supports the 
above mentioned observations that the CH modes seems 
to behave independently of the CC modes. 

7.2. The feature assignments and carriers 

In general, the peak position of the PAH vibrations may be 
affected by charge (anion, cation, and neutral states), size, 
symmetry and molecular structure and heterogeneity. It is 
therefore important to understand the precise transitions 
involved in producing these interstellar features if we are 
to fully interpret what these transitions tell us about the 
nature of the carriers, their history and the local physical 
environment. 
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7.2.1. The 6.2 feature 

The observations presented here (Figs. p| through pi]) show 
that the peak position of the intersteUar PAH CC stretch- 
ing band varies between 6.2 and 6.3 ^m. The intersteUar 
feature seems more red-shaded (i.e. asymmetric with a red 
wing) when it peaks at 6.2 /im (Class A emitters, Table 
2) then when it peaks near 6.25 /im (Class B emitters, 
Table 2). For Class A emitters, the observed half- widths 
(0.11 /im, 28 cm~^) are close to the expected intrinsic line 
width of an emitting PAH molecule or cluster (roughly 30 
cm~^), for Class B emitters the observed widths (0.15 /im, 
38 cm~^) are slightly larger and Class C emitters have the 
smallest widths (0.099 /zm, 25 cm~^). 

The variation in the observed peak positions and pro- 
files are interwoven and may reflect the effects of a class 
of PAHs, each with its individual peak position, or the 
effects of anharmonicity inherent to the emission process 
in hig her vibrationally excited PAHs, or both. 



Anharmonicity 

It is also possible to interpret the observed (red shaded) 
profiles in terms of emission from a single PAH (or a col- 
lection of PAHs with very similar peak position) which is 
(are) highly vibrationally excited. If a molecule is suffi- 
ciently highly vibrationally excited, emission from levels 
above the first excited state become important. Due to 
the anharmonic nature of the potential well, these band 
spacings become smaller and smaller the higher up the vi- 
brational ladder one samples and emission between these 
levels falls slightly and progressively to the red producing 
a long wavelength wing reminiscent of the observed wing 
(Barker et al. 1987). In addition, anharmonic coupling of 
the emitting mode with other modes also shifts the peak 
position of the emitting band to lower energies. Integrat- 
ing over the energy cascade as the emitting species cools 
down will then in a natural way give rise to a red shaded 
profile (Barker ct al. 1987; Pech et al. 2001; Verstraete 



et al. 2001) 



Emission by a PAH family 



As discussed in Sect. <xl and illustrated in Fig. pure 
PAHs (i.e. those containing only carbon and hydrogen) 
composed of > 20-30 C atoms can reproduce the position of 
the interstellar 6.3 /im component but cannot account for 
the 6.2 /im component. Based on the currently-available 
experimental and theoretical data, the most probable car- 
riers of the 6.2 /im component are hetero-atom (e.g. N, 
O, Si) substituted PAHs, PAH clusters and/or PAH com- 
plexes with a metal atom such as iron (i.e. mctallocenes) . 
Interactions in such species may alter a PAH molecule 
in two ways : they may lower the molecule's symmetry 
and they may alter the electron distribution within the 
molecule. Both of these effects may shift the dominant 
CC stretching modes toward shorter wavelengths. 

In summary, the observed profiles of the 6.2 //m band 
complex can be a " blend" of emission bands by different 
carriers with slightly shifted peak position for the main CC 
mode. Pure PAHs represent the 6.3 fim component and 
substituted PAHs or PAH-like species emits the 6.2 /im 
component. The profile of the 6.2 /im band complex re- 
flects then the relative abundance of the different species 
in the PAH population and so different sources contain 
different amounts of each type of species. Typically, the 
PAH population in Post-AGB objects, PNe and isolated 
Herbig AeBe stars is skewed toward pure-C PAHs, while 
HII regions, RNe, non-isolated Herbig AeBe stars and 
galaxies have a dominant contribution from substituted 
complexed PAHs. However, the PAH family in some PNe 
also contains an important substituted PAH population. 
NGC 7027 is a case in point. This interpretation raises 
interesting questions on the origin and evolution of cir- 
cumstellar and interstellar PAHs which will be addressed 



Pech et al. (2001) and [Verstraete et al. (200lD have 
modelled the IR emission spectrum of PAHs based upon 
extensive laboratory studies of the shift in peak position 
as a function of temperature of the emitter which is a di- 
rect measure of the anharmonicity ( Joblin et al. 1995 ). 
They obtained excellent fits to the red shaded appear- 
ance of the observed profiles of the 3.3, 6.2 and 11.2 fj,m 
bands. However, for the 6.2 /im band no good fit to the 
peak position is obtained for sources with class A spec- 
tra, because the laboratory and theoretical studies have 
been limited to pure-C PAHs. Nevertheless, the principle 
remains the same and small, highly vibrationally excited, 
N-substituted PAHs are also expected to have red shaded 
emission profiles. 

Within this interpretation, the class A profiles which 
are highly asymmetric are due to a relatively limited 
number of small highly excited N-substituted complexed 
PAHs. The class C profiles, which peak at 6.3 /im and 
are fairly symmetric (cf. Fig. ^), are carried by much 
less highly excited pure-C PAHs. The "low" excitation 
of these pure-C PAHs may reflect either the cool nature 
of the illuminating source in these two Post-AGB ob- 
jects (CRL 2688 and IRAS 13416) or it may reflect an 
on-average larger size of the emitting species, or both. 
Class B profiles show a less pronounced blue rise and a 
less pronounced red wing. We note that class B sources 
typically have a higher I11.2/I12.7 intensity ratio than 



class A sources (Hony et al. 2001). A high I11.2/I12.7 



in Sect. 7.3. In principle, the family could consist of two 
(distinct) species only. However, we consider that unlikely. 



intensity ratio indica tes the dominanc e of rather large 
150 C-atoms) PAHs ( Hony et al. 2001 ) and anharmonicity 
effects are expec ted to be smaller for such large PAHs 
( iPech et al. 2"00l|) . The YSO, BD+40 4124, which belongs 
to class A, shows a Gaussian profile, possibly indicating 
that the photochemical survival of the fittest (i.e. most 
robust) members of the PAH family has been of great 
importance in this source (Van Kerckhoven et al. 2002, 
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in prep.). 

The 6.2 fj,m band of many of these sources have very 
similar profiles, irrespective of the harshness of the illu- 
minating radiation field - as measured by its strength or 
effective temperature. As exemplified by the model study 
of Verstraete et al. (2001), this indicates that the typical 



size of the emitting PAH is larger in regions which are il- 
luminated by hotter stars. This coupling between size and 
the "colour" of the illuminating radiation field may be 
a natural consequence of emission from a family of PAH 
species; that is, the smallest size of PAHs which can sur- 
vive in a radiation field will depend on the average photon 
energy in the illuminating FUV field. To phrase it differ- 
ently, both the minimum size and the profile of the 6.2 /im 
band may be a measure of the avera ge excitation of the 
emitting PAH ( Verstraete et al. 2001 ) . The more symmet- 
ric profiles of class C may then reflect that the PAHs in 
these cool Post-AGB objects have not yet been exposed to 
harsh radiation fields and their composition still refiects 
the condition during their formation at high temperatures. 
As a corollary, this implies that there are no small PAHs 
(< 25 C atoms) with peak positions longwards of 6.3 /im. 
The coolness of the radiation field and symmetric profiles 
imply then that the PAHs in these sources are only mod- 
erately excited. 

7.2.2. The 7.7 /im complex 

As discussed in the spectroscopy section (Sect. |6.2| ), the 
observed 7.7 nm profiles of class A' can be remarkably 
well reproduced by either pure-C or N substituted PAHs. 
The dominance of the 7.8 /im band in class B' profiles, 
however, remains an enigma. 

When investigating other possible carriers for the 7.8 
/xm component, one should bear in mind the following 
observational facts. First, the 7.7 /zm complex is always 
observed together with the other UIR bands. Second, the 
position of the 7.8 fim component correlates with the ob- 
served intensity ratio ly.e/Ir.s- Third, the strength of the 
7.8 /xm component is correlated with the strength of the 
6.2 fj,m feature. In addition, the different classes of the 6.2 
and 7.7 /im features are directly linked with each other. 
Hence, the carrier of the 7.7 /im complex and in particular 
of the 7.8 /im component should be related to the carrier 
of the other UIR bands. 

Other carriers have been proposed to explain the UIR 



bands. Proposed solid state carriers are QCC (Sakata et al 



1984 ) , soot ( AUamandola et al. 1985|), coal ( Papoular et al 



1989| ), HAG (iColangeh et al. 1995|; ^cott et al 1997|) a^ 



nano-sized carbon grains ( Herlin et al. 1998| ; gchnaiter 



et al. 199£ ). Laboratory measured spectra of the solid 
state materials all resemble the global appearance of the 
observed UIR spectrum. Looking in detail, however, they 
do not match the observed peak positions, the observed 
widths and the observed profiles. In addition, these grains 



would be generally too cool to emit efficiently in the mid- 
IR. In summary, to date the molecular carrier of the so- 
called 7.8 /im component remains unidentified albeit that 
it likely has a highly aromatic character. 



7.2.3. The 8.2 /im feature 

Two Post-AGB stars (GRL 2688 and IRAS 13416) in our 
sample show a peculiar IR spectrum. Instead of a 7.7 fim 
and 8.6 /im feature, they exhibit a broad 8.2 /im feature. 
The other UIR bands are also slightly different. The 3.3 
/im feature has a similar peak position as HD 44179 (type 
2 in Tokunaga et al. 1991). But, it is broader than in 
HD 44179 or in HII regions. In addition, these objects emit 
a symmetric 6.3 /im feature. Unfortunately, the 11.2 /im 
feature is too weak to define its profile and peak position. 
Both sources show a 3.4 /im band. 

There are several ways to interpret these spectra. First, 
the observed spectra could be a combination of emission 
by PAHs, giving rise to a (slightly modified) UIR spec- 
trum, and by an unknown carrier which produce exclu- 
sively the 8.2 /im feature. In this interpretation, the 7.7 
and 8.6 /im bands would be hidden in the strong 8.2 
/im feature. Second, the carriers of the features in these 
sources might have similar GH modes as PAHs but their 
GG modes are different. Since dust is formed in the out- 
ffows of Post-AGB objects, the spectra of these sources 
might then refiect that of freshly synthesised PAHs, dust 
and intermediate compounds. 

This broad 8.2 /im band may well be present in some 
other sources as an underlying plateau (Sect. U). Studies 
of the spatial distribution of this plateau have shown that 
it is carried by a component which is independent of that 



of th e 6.2 and 7.7 /im bands ( Bregman 1989 ; Cohen et al 



1989 ). Energetic arguments suggests that the carrier of 
this plateau contains ~400 G-atoms and, hence, the carrier 



may be in the form of PAH clusters ( Bregman 1989 ). 

In these two Post-AGB objects, the carrier of the 8.2 
/im band may also be in the form of larger grains. The 
dust (and gas) in these two Post-AGB objects is very close 
to the central star and hence dust particles might attain 
high enough temperature to emit strongly around 8 /im. 
The large width of the 8.2 /im band lends some credence 
to a grain-like carrier for this band. Various carbonaceous 
materials show an emission near 8 /im, including HAG, 



QGG, coal, and partially hydrogenated Ggo ( 


Mortera & 


Low 1983; pakata et al. 1984; polangeh et al. 1995; 


Guil- 


lois et al. 1996|; 3cott et al. 1997; 


Schnaiter et al. 1999; 


Stoldt et al. 2001). However, the profiles in these solid 



state materials peak close to 8 /im and are much broader 
then the 8.2 /im component. These materials do provide 
reasonable fits to t he 8 /im feature as observed toward 



IRAS 22272-1-543 5 ( |Buss et al. 1993| ; puillois et al. 1996 



Kwok et al. 2001 ) but cannot be the carriers of the 8.2 /im 
feature. 
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7.3. Formation/Evolution 

It was recognised some time ago that there are two main 
classes of contributors to the 7.7 /im band, one peaking 
near 7.6 /im which is associated with HII regions, and 
one p eaking near 7.8 /im, associated with planetary neb- 
ulae ( Brcgman 1989| ; [Cohen et al. 1989|) . In this paper, 
we reported three main classes. These classes are indeed 
related with the type of object (see Sect. ^). In general, 
HII regions, RNe, non-isolated Herbig AeBe stars and the 
extragalactic sources form one class (A) with dominant 
bands at ^ 6.22 /um, ~ 7.6 /um and ^ 8.6 /im. The local 
radiation field Gq in the latter sources ranges from 3E2 to 
7E6. Isolated Herbig AeBe stars, PNe, HD44179, HR4049 
- with a local radiation field Gq between 6E4 and 2E7 - 
form a second class (B) with dominant bands in the range 
of 6.24-6.28 /im, longwards of 7.7 /imand longwards of 
8.62 /im. Two peculiar Post-AGB stars form a third class 
(C) with a 6.3 /im and a 8.22 /im feature; with a local 
radiation field Go of 5E3 in one of these two stars. Two 
Post-AGB stars and one PN belong as well to class A. 
As discussed in the previous section, class A is proba- 
bly built up by nitrogen substituted or complexed PAHs 
while class B seems to be dominated by pure PAHs in 
their dust collection. The PAH spectrum in the 6 to 9 /im 
region apparently reflects local physical conditions and the 
accumulated effect of processing from the formation sites 
in the AGB or post-AGB phases to the ISM. 

The two main PAH classes identified through the peak 
position of their CC modes - classes A and B - imply 
two distinct histories. Because these two PAH classes are 
connected to classes of astronomical objects, these his- 
tories are likely locally determined. Hence, we interpret 
class B as the PAHs formed in the stellar ejecta pre- 
sumably through chemical processes similar to terrestrial 
soot formation. Extensive laboratory studies and theoret- 
ical calculations have shown that in that case highly con- 
densed, "pure" carbon PAHs are the favoured molecular 



intermediaries in the dust condensation route (Frcnklach 



1990| ; [Frcnklach fc Feigclson 1989| ; |Chcrchncff ct al. 1992|) 



The PAHs in class A, on the other hand, likely repre- 
sent a modification of those in class B in the harsh en- 



vironment of the interstellar medium ( AUamandola ct al 
1999| ). Interstellar PAHs could be substantially chemically 



processed in the warm gas of strong shock waves. Alter- 
natively, energetic processing through UV and/or cosmic 
rays may lead to some modification of the PAH structure 



(3trazulla et al. 1995 


2002a 


b|). 



Hony et al. (2001) concluded from their study of the 
CH out-of-plane bending modes that evolution is an in- 
tegral aspect of the life of interstellar PAHs. Specifically, 
the spectra of PNe point towards the presence of large 
(~ 150 C-atoms) compact PAH with regular molecular 
edge structures. In contrast, HII regions are dominated by 
highly irregular molecular edge structures. In their view. 



these molecular differences are driven by energetic pro- 
cessing in the harsh condition of the ISM of PAHs initially 
formed in stellar outflows. Our study of the CC modes 
complements this view. The spectral variation in the CC 
modes in these very different environments similar attest 
to the chemical evolution of the family of PAHs. 



8. Summary 

The most striking aspect of the features in the 6-9 /im 
region is their variability. All features shift in peak position 
from source to source, show different profiles and each 
seems to be composed of several subfeatures. Moreover, 
the variations in the 6.2 and 7.7 /im bands seem to be 
correlated with each other. In addition, these variations 
depend on the type of source considered and apparently 
reflect local physical conditions or the accumulated effect 
of processing from the formation sites in the AGB or post- 
AGB phases to the ISM. In particular, the sources with 
a profile A 6.2 /im feature have a "7.7" feature peaking 
at 7.6 /im, while for those with a component B 6.2 /im 
feature, the "7.7" feature peaks longwards of 7.7 /im. The 
class C objects, with a 6.2 /im feature peaking at 6.3 
/im do not show a "7.7" feature but instead show a broad 
emission band at 8.2 /im. These variations stand in marked 
contrast to the behaviour of 3.3 and 11.2 /im bands whose 



profil es and peak positions are quite invariable ( Hony ct al 
^00 1| , Peeters et al. 2002, in prep.; van Diedenhoven et al. 
2002, in prep.). 

We have summarised laboratory data and quantum 
chemical calculations for the 6~9 /im region. We attribute 
the observed 6.2 /im profile and peak position to the com- 
bined effect of a PAH family and anharmonicity with pure 
PAHs representing the 6.3 /im component and substi- 
tuted/complexed PAHs the 6.2 /im component. The 7.6 
/im component is well reproduced by both pure and sub- 
stituted/complexed PAHs but the 7.8 /im component re- 
mains an enigma. In addition, the exact identification of 
the 8.22 /im feature remains unknown. 

The observed spectral variations in the CC modes are 
coupled to the astronomical characteristics of the sources 
(object type). We find this strong support for the presence 
of a family of PAHs whose composition and/or emission 
characteristics are sensitive to the local physical condi- 
tions. An analysis of the CH out-of-plane bending modes 
in a similar sample has drawn essentially the same con- 
clusion. Apparently, the interstellar PAH family is readily 
processed in space environments. 

The past decade has witnessed great experimental and 
theoretical progress in understanding the spectroscopic 
properties of PAHs under interstellar conditions. However, 
the new observational data presented here poses signifi- 
cant new questions concerning the nature of the band car- 
riers, questions whose answers will yield insight into the 
nature of the emitters and history of the emission zones. 
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